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ABSTRACT 

Gamma-Ray Bursts (GRBs) are likely to have made a number of significant 
impacts on the Earth during the last billion years. The gamma radiation from 
a burst within a few kpc would quickly deplete much of the Earth's protective 
ozone layer, allowing an increase in solar UVB radiation reaching the surface. 
This radiation is harmful to life, damaging DNA and causing sunburn. In addi- 
tion, N0 2 produced in the atmosphere would cause a decrease in visible sunlight 
reaching the surface and could cause global cooling. Nitric acid rain could stress 
portions of the biosphere, but the increased nitrate deposition could be helpful 
to land plants. 

We have used a two-dimensional atmospheric model to investigate the effects 
on the Earth's atmosphere of GRBs delivering a range of fluences, at various 
latitudes, at the equinoxes and solstices, and at different times of day. We have 
estimated DNA damage levels caused by increased solar UVB radiation, reduction 
in solar visible light due to NO2 opacity, and deposition of nitrates through 
rainout of HNO3. For the "typical" nearest burst in the last billion years, we find 
globally averaged ozone depletion up to 38%. Localized depletion reaches as much 
as 74%. Significant global depletion (at least 10%) persists up to about 7 years 
after the burst. Our results depend strongly on time of year and latitude over 
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which the burst occurs. The impact scales with the total fluence of the GRB at 
the Earth, but is insensitive to the time of day of the burst and its duration (Is to 
1000s). We find DNA damage of up to 16 times the normal annual global average, 
well above lethal levels for simple life forms such as phytoplankton. The greatest 
damage occurs at low to mid latitudes. We find reductions in visible sunlight 
of a few percent, primarily in the polar regions. Nitrate deposition similar to or 
slightly greater than that currently caused by lightning is also observed, lasting 
several years. We discuss how these results support the hypothesis that the Late 
Ordovician mass extinction may have been initiated by a GRB. 

Subject headings: astrobiology - gamma rays: bursts - atmospheres 

1. Introduction 

Knowledge of the effects of gamma-ray bursts (GRBs) on the Earth and life is important 
for establishing where and when life is likely to exist in our own Galaxy and others. Concepts 
such as the Galactic Habitable Zone (Gonzalez et al. 2001) depend on information about the 
rate and distance of effectiveness of GRBs, along with other transient radiation events such 
as supernovae. In addition, GRBs may have played a role in mass extinctions in the Earth's 
history (e.g. Melott et al. (2004) and references therein). 

As discussed in Melott et al. (2004) and Dermer & Holmes (2005) it is likely that 
a GRB has occurred at least once in the last Gy close enough to have dramatic effects on 
stratospheric ozone, leading to detrimental effects on life through increases in solar ultraviolet 
(UV) radiation which is strongly absorbed by ozone. 

In Thomas et al. (2005) we presented results of an initial set of parameters modeling 
the impact of a gamma-ray burst on the earth's atmosphere and subsequent effects on life. 
In this study, we present further results for various parameter values. In particular, we vary 
incident fluence, time of year of the burst, impact latitude, time of day, and duration of the 
burst. 



2. Methods 

2.1. Atmospheric Model 

Our modeling of the atmospheric effects of a GRB was performed using the Goddard 
Space Flight Center two-dimensional atmospheric model, which was also used in the study 
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reported in Thomas et al. (2005). We briefly describe the model here, and give a more 
detailed description in Appendix A. 

The model's two spatial dimensions are altitude and latitude. The latitude range is 
divided into 18 equal bands and extends from pole to pole. The altitude range includes 
58 evenly spaced logarithmic pressure levels (approximately 2 km spacing) from the ground 
to approximately 116 km. A lookup table is used for computation of photolytic source 
term, used in calculations of photodissociation rates of atmospheric constituents by sunlight 
(Jackman et al. 1996). Winds and small scale mixing are included as described in Fleming 
et al. (1999) (see Sect. A.2). 

We have employed two versions of the atmospheric model. One is intended for long 
term runs (many years) and includes all transport mechanisms (e.g., winds and diffusion); 
it has a time step of one day and computes daily averaged constituent values. The second is 
used for short term runs (a few days) and calculates constituent values throughout the day 
and night, but does not include full transport. Previously, this version has been used with 
a time step of 225 seconds (Jackman et al. 2001). In the current study we have used a time 
step of one second in order to allow for inputting our GRB ionization over several time steps 
and to properly account for chemistry during and after the burst (see also Sect. A. 3). 

We have not included the effects of any ultra high energy (> 10 18 eV) cosmic rays 
from a GRB (Dermer & Atoyan 2004; Waxman 2004a,b; Dermer & Holmes 2005), due to 
uncertainty as to whether and at what energies GRBs may produce such particles. 

2.2. GRB Ionization Rate Profiles 

Though there is much variation among GRBs, estimates of the typical isotropic equiva- 
lent power of low redshift GRBs range from 4.4 x 10 44 W (Guetta et al. 2005) to 6 x 10 44 W 
(Lloyd- Ronning et al. 2002). Therefore, for the purposes of this study, we assume a "stan- 
dard" burst with power 5 x 10 44 W (isotropic equivalent). Burst durations (considering only 
"long" GRBs) range from 1.66 to 278 s (Preece et al. 2000). We assume for this study a 
standard burst duration of 10 s. We also investigate bursts of 1, 100 and 1000 s duration. 
We do not consider short bursts due the lack of data on rate, luminosity, distance, etc. 

In Melott et al. (2004) we estimated conservatively that the probable nearest burst 
pointed at the Earth in the last billion years is approximately 2 kpc. For our standard burst 
this corresponds to a fluence of 100 kJ/m 2 . Dermer & Holmes (2005) estimated this distance 
to be 1 kpc. Our estimate of the nearest probable burst distance in the last billion years is 
arrived at by scaling from the universal observed rate (as of 2003), taking into account the 
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universal star formation rate, and using the blue luminosity surface density of the Galaxy. 
We then use the fluence from our standard burst at this distance to compute atmospheric 
and other effects. 

Ionization rate profiles due to a gamma-ray burst are input to the model as production 
sources of NO y and HO x . It is assumed that for each ion-electron pair produced, 1.25 NO y 
molecules are produced at all pressure levels (Porter et al. 1976) and 2.0 molecules of HO x 
are produced below 75 km and less than 2.0 (from 1.99 to 0.0) for altitudes greater than 75 
km (Solomon 1981). 

The ionization rate profiles are calculated separately from the atmospheric model, fol- 
lowing the method used in Gehrels et al. (2003). This section draws on that publication 
as well as an unpublished manuscript by Claude Laird. In the Gehrels et al. (2003) study, 
gamma-rays were included using the spectrum of SN 1987A. In the present study, the gamma- 
ray differential photon count spectrum used is that of Band et al. (1993), which consists of 
two smoothly connected power laws: 

m( E ) = A (E_y (-E 
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where we take typical values E = 187.5 keV, a = —0.8, (3 = —2.3 (Preece et al. 2000). 
We have not taken variability of these quantities during the burst into account. The flux 
constant, A, is used to scale the total received energy to our desired value (in this study, 
corresponding to fluences of 10 kJ/m 2 , 100 kJ/m 2 and 1 MJ/m 2 ). 

The total photon flux in each of 66 evenly spaced logarithmic energy bins, ranging 
0.001 < E < 10 MeV, is obtained by integrating the above spectrum for each bin. The high 
and low energy regions of the spectrum are integrated separately; the E > (a — j3)E region 
is integrated analytically and the E < (a — (3)E energy region is numerically integrated 
using the "qsimp" routine from Numerical Recipes (Press et al. 1992). 

The total incident energy flux in a monoenergetic beam at the top of the atmosphere is 
Fj° = Ni°{Ei), where iV; is the incident photon flux (calculated from the Band spectrum as 
described above). The photon flux is propagated vertically through a standard atmosphere 
(adjusted for the appropriate latitude and time when input to the atmospheric model), and is 
attenuated with altitude by an exponential decay law, with frequency-dependent absorption 
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coefficients taken from a lookup table (Plechaty et al. 1981). 

The photon flux for the ith energy bin which remains at the jth altitude layer is given 
by Nij = Ni°e~^ iXj , where Xj is the column density (in g/cm 2 measured from the top of the 
atmosphere and //; is the absorption coefficient (from Plechaty et al. (1981)). The photon 
flux deposited in the jth layer with energy (Ei) is AJVjj = Nij-i — N iy j. The corresponding 
energy flux deposited is F it j = AN i j(E i ) (in MeV/cm 2 /s). The total ionization rate is the 
sum over all energies: 

1 66 F ■ 

qtot > 3 = 35eV^£% (2) 

i=l J 

where 35 eV is the energy needed to produce one ion pair (Porter et al. 1976) and AZj is 
the thickness of the jth altitude layer. 

The vertical ionization rate profiles are mapped onto an altitude-latitude grid for use 
by the atmospheric model as follows. First, a series of ionization rate profiles as functions 
of altitude and latitude are calculated for a flat Earth for incidence angles of 5° - 85° 
from zenith, corresponding to latitude bands 10° wide centered on 85° - 5°, respectively, 
for the hemisphere illuminated by the GRB at zenith point 90° N. These profiles are then 
interpolated and the zenith frame latitudes are transformed into standard Earth coordinates 
(longitude and latitude) according to the impact latitude chosen for the event. Finally, a 
zonally averaged profile (as a function of altitude and latitude) is produced by averaging 
over all longitudes for each 10° latitude band. This profile is then input to the atmospheric 
model as a production source of NO y and HO x , as described above. 



2.3. Simulation Runs 



Our simulation runs are performed as follows. Initial conditions are obtained from long- 
term runs (roughly 40 years) intended to bring the model to equilibrium. These runs end 
at the various times of year at which we input the GRB. Constituent values from these 
equilibrium runs are read in by the 1 s time step version of the model which runs for 7 days, 
beginning and ending at noon. Ionization due to the gamma-rays (as calculated off-line, see 
Sect. 2.2) is input as a step function (top hat) with duration 10 s. Most runs input the 
burst at noon on day 4. Bursts at other times of day occur after this. We use day 4 as the 
input point in order to allow the short-term model to come to equilibrium. This is necessary 
as the two versions have different time scales for chemistry (1 s versus daily averages) and 
differences in transport. The run then continues for three more days to allow for good time 
resolution of the chemistry directly following the burst. Running for much longer than 7 
days would introduce errors due to the limited transport in the short term model. Hence, at 
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the end of 7 days we then feed constituent values to the long term (1 day time step) version 
of the model which runs for 20 years, by which time the atmosphere has recovered from the 
GRB perturbation. Runs are also performed using the same steps without input of GRB 
ionization, in order to make comparisons between perturbed and unperturbed runs (e.g. 3 
percent change). 

2.4. Biological Effects 

In order to quantify some of the effects on life due to a GRB impact, we have computed 
the direct effect of the increased solar UVB (due to depletion of ozone) by determining the 
irradiance in this waveband at the surface of the Earth and convolving it with a biological 
weighting function which describes the effectiveness of particular wavelengths in damaging 
DNA molecules. 

In order to determine the solar UVB reaching the surface, we take the irradiance at 
the top of the atmosphere (obtained from data taken by the SUSIM ATLAS instrument, 
available: http : //wwwsolar.nrl. navy.mil/susim. html) and determine the irradiance at the 
surface (for a given wavelength band) using the Beer-Lambert Law, 

I = lQe -a(X)N/cose (3) 

where <r(A) is the absorption cross section of ozone as a function of wavelength, A, (e.g. 
Yoshino et al. (1988)), N is the vertical column density of ozone (taken from our modeling 
results) and 9 is the solar zenith angle (calculated for a particular latitude, day of year and 
time of day) (Madronich 1993). Only ozone absorption effects are included in this calculation 
since the effect of scattering at these wavelengths is comparatively small (see Sect. 4.0.7). 
The surface irradiance is then multiplied by the value of the biological weighting function 
(Setlow (1974); Smith et al. (1980); see Fig. 1) for that wavelength band and this combined 
DNA damage value is then integrated over the entire wavelength range. This value is then 
scaled by cos# to account for insolation. The above calculation is done for each hour and 
then averaged over 24 to get a daily average DNA damage value for a given latitude and 
time point. 

3. Results 

The primary atmospheric results of our simulations are increases in NO y and decreases 
in O3. Ozone column densities computed in the model are then used to calculate the UVB 
flux at the Earth's surface and the resulting DNA damage. For purposes of reference and 
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validation of our model's ability to correctly simulate normal ozone variations we show in 
Fig. 2 a plot of 3 column density for the unperturbed atmosphere as a function of latitude 
and time, over one year, starting in January Units shown are both Dobson units and 
10 18 cm~ 2 . A Dobson unit describes the thickness of a column of ozone at STP and is 
defined as 1 DU = 0.01 mm thickness (or, 1 DU = 2.69 x 10 18 cm -2 ). CFC's and other 
anthropogenic constituents are not included in our runs (see Sect. A), so the values and 
distribution of 3 in this plot differ somewhat from those in the present day atmosphere. 
Anthropogenic constituents are excluded because we are primarily interested in applying our 
results to pre-industrial time periods. 

Results for an initial set of parameters were reported in Thomas et al. (2005). We repeat 
those results here, in the context of a wider exploration of parameters. First, we discuss the 
effects of varying the time of year and the latitude at which the burst occurs for a given 
fluence (Sect. 3.2). Then, we discuss the effects of the burst occurring at different times of 
day and with various durations, for a given fluence (Sect. 3.3). Next, we discuss the effects 
of varying the burst fluence (Sect. 3.4). We also present results quantifying the impact of 
reduced O3 levels on life (Sect. 3.5). Finally, we discuss the reduction of sunlight due to 
increased NO2 levels and the deposition of nitrates due to rainout of HNO3 (Sect. 3.6). Note 
that most figures presented here are available in color in the online version of the paper. 

3.1. The Importance of Sunlight 

As will be seen throughout the discussion that follows, the presence or absence of sunlight 
at a given location and time has important implications for our atmospheric results. This 
is due to photolytic reactions involving 3 and constituents of NO y , most importantly NO 
and N0 2 . 

The primary creation mechanism for O3 is by dissociation of O2 by solar UV at wave- 
lengths below 242 nm. Higher wavelengths dissociate O3, which is why ozone provides an 
effective shield at these wavelengths. A consequence of the photo-production of O3 is that 
ozone is primarily created in the tropics, where sunlight is most constant and intense. Of 
course, transport in the upper atmosphere is primarily from the tropics to the poles, and so 
ozone is moved to the poles after being produced near the equator. Longitudinal transport 
is quite efficient and only takes a few weeks to circle the globe. 

Photolysis is also important for NO y constituents and their effects on ozone. NO and 
NO2, the primary ozone depleting components of NO y , are dissociated by sunlight at wave- 
lengths below 191 and 400 nm, respectively. In addition, other NO y constituents are pho- 
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tolyzed to produce NO and N0 2 . (See Table 11 for photolysis reactions and Table 12 for 
solar fluxes included in the model.) Photolysis of NO has two important implications for 
both NO y and O3. This reaction yields atomic N which can reduce NO through the reac- 
tion N + NO — > N 2 + O, thereby decreasing concentrations of NO y . In addition, atomic O 
is produced which can actually increase O3 through reaction with 2 . 

The diurnal cycle of NO and N0 2 exhibits an interchange between these compounds. 
NO concentrations decrease as the sun sets due to the removal of its photolysis source and 
its rapid oxidation to N0 2 . During the day, NO is created partly through photolysis of NO2 
thereby reducing the amount of NO2 present. NO2 is more strongly affected by photolysis 
than NO, since its upper wavelength limit (400 nm) is higher than that for NO (191 nm) 
and there are more photons present in the solar flux below 400 nm than 191 nm (Table 12). 

These interchanges can be seen in Fig. 3 where we plot column densities of several 
constituents at —5° latitude, normalized to their maximum values over this time period. 
This data is for a 100 kJ/m 2 burst over the equator in late March at noon. Time indicates 
the burst input. It is apparent from this plot that the burst causes an immediate increase in 
the concentration of NO y (primarily through NO). After the burst, normal diurnal processes 
dominate the relative amounts of these compounds at any given time, though at much higher 
levels than normal, and with the addition of the long term depletion of O3. 

The effects of sunlight on NO, N0 2 and 3 are important for understanding many 
features of our results and we will refer to them throughout the following discussion. 

3.2. Varying Time of Year and Latitude 

We will now discuss results similar to those presented in Thomas et al. (2005) for bursts 
occurring at the equinoxes and solstices, over five different latitudes. Figures 4, 5, 6, and 
7 show results for bursts of fluence 100 kJ/m 2 input in late March (repeating the results 
from Thomas et al. (2005)), June, September and December, at latitudes +90°, +45°, 0° 
(equator), —45°, and —90°. For all these cases the burst is input at noon. In each plot, the 
burst is input at time 0. 

Figure 4 shows the vertical column density of NO y , in units of 10 16 cm -2 , at each 
latitude point in the model over time for each combination of event time of year and latitude. 
Similarly, Fig. 5 shows the vertical column density of O3. Included are scales in both Dobson 
units and 10 18 cm -2 . Figure 6 shows the percent difference at a given location and time 
(between a run with gamma-ray input and one without) in vertical column density of O3. 
(Note that in this plot white indicates values equal to or greater than 0.0. The "wedges" 
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of white that appear in the polar regions in some frames are the only areas of O3 increase. 
This will be discussed further below.) The full long-term duration of the effects can be seen 
in Fig. 7, which shows the global average percent change in ozone as a function of time for 
15 years after the burst (with one year of change before the burst). We summarize the 
maximum depletion of ozone (globally averaged and at a given latitude) in Tables 1 and 2. 

To give a sense of the distribution of atmospheric effects with altitude we also include 
Fig. 8 and 9, which show pointwise percent difference in NO y and O3 between a run with 
a burst and one without, at each altitude and latitude at one year after the burst occurs. 
These plots are for a burst in late March over the equator. 

An interesting feature apparent in Fig. 9 is that while ozone levels show a decrease at 
altitudes above about 20 km, they are actually increased at or below this altitude, particularly 
around the equator. (The white "boundary" starting at about 15 km in the polar regions and 
moving upward toward the equator separates areas of O3 decrease at higher altitudes and 
O3 increase at lower altitudes.) This production at lower altitudes is due to the following. 
Normally, ozone is concentrated between 30 and 40 km altitude and is sparse at lower 
altitudes. However, the burst depletes ozone primarily in this altitude range, which allows 
solar UV (normally absorbed by the ozone) to penetrate to lower than normal altitudes. 
This UV then creates ozone (by photolysis of O2) at altitudes where normally there is little 
or none. (Jackman & McPeters (1985) discuss observation of this effect in the case of a solar 
proton event.) This effect is particularly strong around the equator, since the intensity of 
solar UV is greatest there. Therefore, since this plot shows a comparison between a run 
without a burst (where little ozone exists at lower altitudes) and one with a burst (where 
UV penetrates more deeply and creates lower altitude ozone), we see an increase at these 
lower altitudes. This effect actually somewhat mitigates the total column depletion. 

Several important features are apparent in Figs. 4-7. First, immediate atmospheric 
effects are seen, i.e. production of NO y and depletion of O3, beginning around the latitude 
over which the event occurs and spreading quickly to other latitudes. Note that for polar 
bursts, even in the long term, effects are isolated to the respective hemisphere. Similarly, 
even for bursts at ±45° the opposite hemisphere experiences much less intense effects. This is 
due to the fact that transport in the atmosphere is primarily poleward from the equator and 
transport across the equatorial line is minimal. Therefore, NO y produced in one hemisphere 
is for the most part not transported to the opposite hemisphere. 

Isolation of effects to a hemisphere means that polar and ±45° bursts result in smaller 
globally averaged ozone depletions as compared to bursts over the equator in most cases (see 
Table 1). This is true even though maximum localized depletions are typically larger for 
polar bursts (Table 2). Local intensity is understandable in these cases, since most of the 



-10- 



NO y is produced directly in the polar regions and subsequently is disbursed over a smaller 
area. This high concentration contained within the polar region allows a greater effect on 
ozone. 

For all cases, ozone depletion is long-lived (Fig. 7). Globally averaged depletions of 
at least 10% are present up to 5 - 7 years and full recovery is not achieved until 10 - 12 
years after the burst. Recovery occurs at the equator first. This is normal in the terrestrial 
atmosphere, as ozone is primarily produced in tropical latitudes and transported poleward. 
Also, as discussed above, NO y produced at the equator is moved toward the poles fairly 
rapidly, so concentrations are significantly reduced at the equator (reducing the depletion of 
ozone there) within a year or so, even for equatorial bursts. 

It is interesting to note that for all the cases discussed so far recovery of ozone takes 
approximately a decade. Why this time scale and not some other? The primary controlling 
factor here is the lifetime of the ozone depleting NO y compounds which are generated by 
the burst. These compounds are removed from the atmosphere mainly through transport 
down to the troposphere where they are rained out. Fleming et al. (1999) and Fleming 
et al. (2001) compare results of the GSFC model with measurements of the "age of air" 
as determined from measurements of long-lived gases such as SF 6 and C0 2 . This "age of 
air" at the tropopause is on the order of 5 years. Only after concentrations of NO y in 
the stratosphere are significantly reduced by this transport process can O3 levels recover, 
a process which is also dominated by transport and hence occurs on a similar timescale to 
the removal of NO y . Therefore, the combination of removal of NO y and replacement of O3 
results in a total recovery time of approximately a decade. 

Greater depletion of ozone in the polar regions is a general feature of these results and 
is observed in the present day atmosphere as well. Figure 6 gives a somewhat exaggerated 
sense of the effect of depletion at the poles because ozone is initially high there (see Figs. 2 5). 
Larger ozone depletions at the poles are due to the preferential transport of NO y poleward, 
combined with the long lifetime of the enhanced NO y in the polar stratosphere. This long 
lifetime is primarily due to lower photolysis rates at high latitudes (from low solar incidence 
angle and dark polar winters; see Sect. 3.1). 

The enhanced NO y in the polar regions, including HNO3, leads to an enhancement of 
nitric acid trihydrate (NAT) polar stratospheric clouds (PSCs). These NAT PSCs facili- 
tate heterogeneous reactions that result in greater ozone depletion by halogen (chlorine and 
bromine) constituents. The contribution of PSCs is especially strong in the south polar re- 
gion where a stronger polar vortex with colder stratospheric temperatures is in place during 
the winter. This leads to the south polar ozone "hole" which appears annually in the present 
day atmosphere. 
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An effect of PSCs is seen in Fig. 4 where column density of NO y is reduced during polar 
winter, especially at the south pole. (The feature is present in the north, but is less intense.) 
This is due to the fact that during polar night most of the NO y is in the form of HNO3 
and, particularly at the south pole, this HNO3 is largely in solid form in PSCs. This form of 
HNO3 is not included in the summation of total NO y and so the NO y column density shown 
here appears reduced during these times. Note, however, that it returns to high values at 
the end of polar night, when solid HN0 3 is converted back to N0 2 . Over the long term, the 
column density of NO y which "reappears" after polar night is less than it was before. This 
is because some PSC particles fall down to the troposphere where HNO3 is rained out. This 
is the primary mechanism whereby NO y is permanently removed from the atmosphere. 

While the particular atmospheric conditions which lead to the north-south asymmetry 
seen here are peculiar to the present-day configuration of continents, the effect of varying 
burst latitude and time of year are much larger. For example, the asymmetry is reversed for 
an equatorial burst in September, and the strength of the asymmetry is less for equatorial 
bursts in June and December. 

The relation of the time of year of the burst to the season at the polar regions is 
an essential factor in the variation of several features between the different burst cases. 
This includes differences such as the shift in north-south asymmetry as well as variation 
of maximum ozone depletion (both globally averaged and localized) for bursts over a given 
latitude which occur at different times of year. Other features are affected by this timing as 
well and will be discussed later in this section. The key variable is how much time elapses 
between when the burst occurs and when polar night ends. Figure 10 shows (in arbitrary 
units) the intensity of sunlight (at noon) at each latitude over 4 years, starting in late March, 
June, September and December, similar to Figs. 4-6. This figure may allow the reader to 
more easily relate the solar intensity at a given time of year and latitude and the effects to 
be discussed below. 

Relative timing of the burst and polar spring is important due to the fact that photolysis 
is an important factor in the chemistry of NO y and ozone and hence the presence (and 
incidence angle) of sunlight is critical (see Sect. 3.1). NO y constituents, particularly NO and 
NO2 (the most important for O3 depletion), are strongly photolyzed. Therefore whether the 
Sun is up in the polar regions, where NO y tends to concentrate, is important for how much of 
an effect on ozone these compounds have. Relative timing of the production of NO y by the 
burst and the appearance of the Sun after polar winter is the primary cause of differences in 
distribution and intensity of ozone depletion for bursts at different times of year. Important 
here is the fact that it takes some time for NO y to be transported to the polar regions when 
the burst is equatorial. Since the removal time for NO y is on the order of several years, these 
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compounds build up in the polar regions. 

For instance, for an equatorial burst in March, high southern latitudes are in darkness 
for 5-6 months after the burst (depending on latitude), whereas at high northern latitudes 
the Sun is up for most of the day during that time period and the next polar spring does not 
begin for some 11-12 months. This means that NO y produced at the equator and transported 
to the northern polar regions may be reduced through photolysis and subsequent reactions 
while that transported to the southern polar regions is not. Hence, the concentration of 
NO y is higher for longer in the high southern latitudes (see Fig. 4). Note that the opposite 
scenario is true for an equatorial burst in September, when the relative timing of the burst 
and spring at the poles is reversed. This may also be seen in the ozone depletion resulting 
from the concentration of NO y (Figs. 5 and 6). 

For bursts which occur in June and December, NO y is somewhat more evenly distributed 
with latitude and ozone depletion effects are reduced, as compared to bursts in March and 
September. For these cases, the time between when the burst occurs and the soonest polar 
spring is smaller for both poles. If polar spring occurs too soon after the burst, significant 
amounts of NO y will not be transported into the region before night returns and so the 
reduction of concentration by photolysis will be lessened. This is apparent for an equatorial 
burst in June at high southern latitudes and in December at high northern latitudes. 

For bursts over the poles and ±45° latitude there are likewise variations in ozone deple- 
tion with burst timing, though primarily in intensity rather than distribution, since for these 
bursts most of the NO y produced is contained within a hemisphere, as discussed earlier. 
Consider a north polar burst. The smallest depletion of ozone occurs for a burst in March. 
This is due to the fact that the Sun is up at most high latitudes at this time and is for several 
months after the burst. For a June burst the depletion becomes somewhat greater, since 
north polar night begins sooner after the burst. A December burst has yet higher depletion; 
here the burst occurs in the middle of north polar night. The greatest depletion occurs for 
a burst in September. In this case, the burst occurs just as north polar night is beginning. 
Note that the relative levels of depletion are reversed between March and September and 
June and December for a south polar burst, simply due again to the relation of the time of 
year of the burst and the season at high southern latitudes. 

Similar relationships are present for the bursts over ±45° latitude. Here there is the 
additional feature of how much depletion occurs in the opposite hemisphere. For instance, a 
March burst at —45° exhibits very little ozone depletion in the northern hemisphere, while a 
September burst results in a fair amount of depletion in the northern hemisphere, but less in 
the southern, as compared to the March case. The same phenomena described above apply 
here as well. In these cases, of course, most NO y is produced and remains in the hemisphere 
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over which the burst occurs, as in the polar cases. However, since the burst occurs at mid 
latitudes, some NO y is produced at the equator or even in the opposite hemisphere. 

A feature which is obvious in Fig. 6 is the presence of brief "spikes" of ozone production. 
These spikes appear with varying intensity for most burst cases. Some occur as lower deple- 
tion than their surroundings, while some exhibit an actual increase in ozone as compared to 
conditions without a burst. All occur close to the poles and are short lived (1-2 months). 
The presence and intensity of these spikes are related to the relative time of year of the burst 
and the season at the polar regions, similar to the effects discussed above. 

Ozone production at the start of polar spring in the presence of high concentrations of 
NO y can be understood as follows. As discussed above, NO y exists primarily as HNO3 during 
polar night. When the Sun suddenly appears, NO2 is produced by photolysis of HNO3. NO2 
is then photolyzed to NO and O and NO is further photolyzed producing additional O. Direct 
production of O3 can then occur when O reacts with O2. In addition, decreasing amounts of 
PSCs (as the temperature rises) and increasing amounts of NO2 reduces the ozone depleting 
action of chlorine and bromine compounds. (N0 2 reacts with CIO and BrO, sequestering CI 
and Br in C10N0 2 and BrON0 2 which do not directly interact with ozone.) 

3.3. Varying GRB Duration and Time of Day 

We will now discuss the differences in results due to varying the time of day at which the 
burst occurs. All the results presented so far are for a burst input at noon. The short duration 
of the burst combined with the normal variation of 3 and NO y constituents (especially NO 
and N0 2 ) over the day-night cycle make it reasonable to expect that bursts which occur at 
different times of day, especially during daylight versus dark, will lead to differing results. 

In Table 3 we present the maximum ozone depletions, globally averaged and localized, 
for bursts occurring every two hours during the day and night, starting at noon, in late 
March, over the equator. As can be seen, there is a difference of 3 to 4 percent between 
daytime and nighttime bursts. However, those which all occur during daylight hours or 
all occur during dark hours are nearly identical to each other. We have also investigated 
a burst at midnight over the equator in late June, which is the maximally different case 
from March, in terms of season. In this case, the maximum globally averaged and localized 
ozone depletion is -40% and -59%, respectively. This compares to -37% and -55% for a 
noon burst at this time of year and latitude. Again, the difference between day and night 
bursts is 3 to 4 percent. The greater ozone depletions for a nighttime burst as compared to 
a daytime burst can be explained by the fact that during darkness NO y can be somewhat 
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more easily produced by a burst and accumulated. (See the discussion above on the effects of 
sunlight on NO y compounds.) For example, the largest column density value for the March, 
equator, burst at noon is 74.23 x 10 16 cm -2 , while for the similar case at midnight the value 
is 74.52 x 10 16 cm~ 2 . 

The long term results for a midnight burst (i.e. results like those in Figs. 4-6) are nearly 
identical to those for a noon burst and so we will not reproduce them here. We do present, 
however, in Figs. 11 and 12 the results of the short term model for NO y and O3 column 
density, showing the immediate effect of the burst. These figures are for a burst at noon, 
in late March, over the equator. The run starts at noon and the burst is input 72 hours 
later (at time in the plots). Note that the scales here are different than in Figs. 4 and 
5, in order to give a clearer sense of the immediate production/destruction, particularly for 
ozone, where maximal effect is not seen for some months after the burst. 

We have also investigated bursts of various durations with total energy fixed. We find 
for bursts of durations between 1 s and 1000 s there is no significant difference in long term 
results. The implication that total energy is more important than duration is in concordance 
with the fact that ozone depletion results for a SN (Gehrels et al. 2003) of the same fluence 
are not greatly different from our results. In that study, a fluence of approximately 15 kJ/m 2 
input over 300 days resulted in maximum globally averaged ozone depletion of about 27%. 
For a GRB delivering the same fluence over 10 s we find a maximum globally averaged ozone 
depletion of 16%. Of course, such a GRB could be 500 times farther away and produce this 
similar effect. 

We would like to understand why total energy input is more important than duration of 
that input. The timescale of our burst range (1-1000 s) is long compared to the ionization and 
dissociation timescales, which are nearly instantaneous. Additionally, the input timescale is 
short compared to the ozone depletion timescale, which is hours to days. To the depletion 
process, all bursts in this duration range are effectively impulsive and provide equal amounts 
of initiating compounds. 

To first order, for the energy regime in which we are working, every photon interacts 
high in the atmosphere by gross ionization and dissocation through many encounters. The 
effect of one photon is essentially independent of another and is proportional to how many 
ionizations and dissociations it can cause, which just depends on its total energy. So, the 
overall effect is, to first order, a result of the total energy of the burst, not its duration 
(within the range we consider here). It is possible that the individual photon "hardness" 
(i.e. the spectrum) may have an effect if it is extremely different because of the atmospheric 
cross section differences. This would, however, be a rather unusual GRB. We are currently 
working on this aspect, as well as including effects from the GRB afterglow. 
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3.4. Varying Burst Fluence 

There is a wide range of possible intrinsic energies and distances from the Earth for a 
burst in the last billion years (see Sect. 2.2 and Melott et al. (2004)). Our primary event 
was selected as a conservative estimate of the probable largest burst in the last billion years. 
We have also modeled bursts delivering two other fluence values to the Earth. 

We will first discuss the results for a burst of incident fluence 10 kJ/m 2 . This fluence 
corresponds to our standard burst at a distance of approximately 6.5 kpc. We will not 
present as complete a set of results as that discussed above for the 100 kJ/m 2 case, since the 
main qualitative features due to time of year, latitude and time of day are the same, though 
the intensity in all cases for this lower fluence is, of course, reduced. We present results for 
a burst in late March over the equator, for comparison with results discussed in Thomas 
et al. (2005) for the 100 kJ/m 2 case. We also present results for a burst in late September 
over +45° latitude. For the 100 kJ/m 2 case, the highest ozone depletions occur for a burst 
in September over the north pole. However, a burst over +45° is more likely than a polar 
burst (due to surface area considerations), so we have chosen this combination as a probable 
upper limit of effects for this fluence case. 

Atmospheric results for a burst of incident fluence 10 kJ/m 2 in late March over the 
equator, and in late September over +45° latitude are shown in Figs. 13 - 16, which are 
similar to Figs. 4 - 7 for the 100 kJ/m 2 cases. Like the previous fluence cases discussed, 
immediate atmospheric effects are seen in production of NO y and destruction of O3. We 
see here again that effects from a burst over +45° latitude are mostly contained within that 
hemisphere. Likewise, greater ozone depletion is seen at the polar regions for both burst 
latitudes and primarily in the southern polar regions for the equatorial burst, as seen for 
the previous cases as well and described in more detail above. Note here that no ozone 
production "spikes" are present, though the depletion is smaller in the same time- frame as 
the spikes in the previous cases. This difference is simply due to the much smaller quantity 
of NO y produced in this case (maximum column densities more than an order of magnitude 
smaller). Again similarly to the 100 kJ/m 2 cases, localized ozone depletion is greater and 
globally averaged depletion is lower for the +45° burst latitude, as compared to an equatorial 
burst. See Table 4 for values of maximum localized and globally averaged ozone depletion. 
In comparison to the 100 kJ/m 2 March, Equator case, this fluence yields only 0.44 times the 
globally averaged ozone depletion and 0.51 times the localized depletion. Similarly, for the 
September, +45° burst this fluence yields only 0.39 of the globally averaged ozone depletion 
and 0.52 of the localized depletion. It is interesting to note that the effects as measured with 
these quantities are only about half as great for a fluence which is an order of magnitude 
smaller. This is likely due to the fact that there tends to be a "freeze-out" of NO production 
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which will limit the relative impact of higher fruence bursts. 

Figure 16 shows the globally averaged ozone depletion over time for both 10 kJ/m 2 cases 
presented. Note that the depletion for the September, +45° case is longer lived (as can also 
be seen in Fig. 15). Significant depletion (greater than 10%) is evident only as long as 2 
years after the burst (1 year for the March, Equator case) and full recovery is achieved by 
about 8 years after the burst in both cases. 

Although our model is probably not capable of correctly handling the chemistry, heating 
and transport for a burst of much higher fluence than 100 kJ/m 2 , we have also modeled a 
1 MJ/m 2 . (See Sects. 4.0.4 and 4.0.6 for a discussion of uncertainties relevant here.) This 
fluence corresponds to a burst at a distance of about 600 pc, which is the closest probable 
burst in the last billion years, based on less conservative assumptions, as discussed in Melott 
et al. (2004). Due to the high level of uncertainty in modeling this fluence, we will not 
present full results. However, for comparison with the above results, our modeling indicates 
that a burst delivering this fluence in late March over the equator yields globally averaged 
ozone depletion of -65% (1.8 times that for the 100 kJ/m 2 case) and localized depletion of 
-80% (1.5 times the 100 kJ/m 2 case). It is again interesting to note that roughly a factor of 
two separates these depletion values for fluences an order of magnitude different. 

It is interesting to note the trends of NO y production and 3 depletion for the three 
fluence cases presented, with the strong caveats of the uncertainty of our modeling of the 
1 MJ/m 2 fluence and the fact that we have only three data points. That said, we find that 
NO y production scales linearly with fluence while the percent change in globally averaged 
column density of O3 scales with fluence as a power law with index ~ 0.3 over this fluence 
range. The slower increase in O3 depletion with fluence as compared with increase in NO y 
production shows that the efficiency of depleting ozone decreases with increased NO y pro- 
duction. This is partly due to increased interaction of NO y with other species, which limits 
the effectiveness both of NO y on ozone as well as other species (such as CI and Br) which 
deplete ozone. There is also higher production of O3 in cases where higher levels of NO y are 
seen. This effect is seen both in total column density (at polar spring, see Sect. 3.2) and 
at altitudes below the normal peak of ozone concentration (see Fig. 9). Additionally, there 
is only so much ozone present to deplete and complete depletion is never reached, but only 
approached asymptotically. For instance, at 32 km altitude, the pointwise percent change 
in ozone one month after the burst is -64% for a 10 kJ/m 2 , -91% for a 100 kJ/m 2 burst 
and -98% for a 1 MJ/m 2 . This shows that the depletion simply cannot scale linearly with 
fluence, even though the NO y production does, since even the middle fluence value depletes 
nearly all of the ozone at that altitude. 
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3.5. Biological Effects 

Ozone normally absorbs about 90% of solar UVB (wavelengths between 290-315 nm) 
radiation. This wavelength range is particularly damaging to organisms since DNA absorbs 
strongly here. Depletion of ozone due to a burst, as we have modeled, will lead to increased 
solar UVB reaching the surface and penetrating a few to 10's of meters in water. As described 
in Sect. 2.4, we have computed the direct effect of this increased UVB on life by determining 
the irradiance at the surface and convolving this with a biological weighting function (Set- 
low 1974; Smith et al. 1980) which describes the effectiveness of particular wavelengths in 
damaging DNA molecules (see Fig. 1). These results are presented in Figs. 17 and 18. Other 
weighting functions also exist, quantifying effects such as erythema (sun burn, e.g. McKin- 
lay & Diffey (1987)), inhibition of photosynthesis (e.g. Cullen et al. (1992)), melanoma (e.g. 
Setlow et al. (1993)), mortality (e.g. Kouwenberg et al. (1999)), etC. While these all vary 
from each other somewhat, the variation is small. Figure 7.1 of Jagger (1985) shows several 
different weighting functions, all of which have similar shapes over our wavelength range. 

Figure 17 shows DNA damage for all burst latitude and time of year cases discussed for 
100 kJ/m 2 fluence (Sect. 3.2). This includes the burst in late March over the equator, which 
was first presented in Thomas et al. (2005). Figure 18 shows DNA damage for the cases 
discussed for 10 kJ/m 2 fluence (Sect. 3.4). We have normalized the plots by dividing the 
damage by the annual global average damage in the absence of a GRB. Note that the scales 
in the figures are not the same, in order to allow for greater visibility of features. Also, white 
areas indicate relative DNA damage values between 1.0 and 0.0 (see color figures online). 
Experiments suggest significant mortality for marine microorganisms should certainly exist 
in areas where this measure exceeds 2. 

Several interesting features may be noted. For many cases, maximum DNA damage is 
evident immediately around the equator and is fairly short-lived. However, for some cases 
(particularly polar bursts), the maximum values occur somewhat later at mid latitudes 
(e.g. around ±30°). In the long term (more than a few months), greatest DNA damage is 
apparent at mid to low latitudes for all burst cases. This heightened damage lasts several 
months at a time, recurring annually. This latitude dependence and annual recurrence is 
due to combination of the 3 depletion effects with the Sun incidence angle, length of day, 
etc. Note that the highest DNA damage values occur for cases with the highest localized 
ozone depletion. 

Like the atmospheric results, the lower fluence DNA damage results show similar qual- 
itative properties to that for the higher fluence (greatest damage initially at the equator, 
long-term enhanced damage at mid latitudes, annually recurring), but less intense. The 
maximum DNA damage for the 10 kJ/m 2 cases occurs with a March burst over the equator 
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and is 0.39 of the maximum for the 100 kJ/m 2 cases (including a March, equator burst). 
Note that this is a smaller fraction than that for ozone depletion (local or globally aver- 
aged) when comparing between fluences for this burst case. This is due to the non-linear 
relationship between DNA damage and ozone depletion (Madronich et al. 1998). Long-term 
enhancement of damage is relatively small and shorter lived here, indicating that damage 
is probably only significant for the first 2 months or so. These features are understandable 
given the relatively small ozone depletions at mid latitudes for this fluence. 

It is important to note that greater DNA damage is evident at mid to low latitudes, in 
the long term, for all cases. One might think that this damage would be countered by a 
greater evolved UVB-resistance in organisms at low latitudes. However, at least for modern 
organisms, there is no evidence that temperate zone phytoplankton are any more UVB- 
resistant than Antarctic plankton (Prezelin 2004). Thus, one might predict that greater 
ecological damage and extinction would be likely near the equator. It is interesting that the 
late Ordovician mass extinction seems to be alone in having recovering fauna preferentially 
derived from high-latitude survivors (Jablonski 2004). 

3.6. Opacity and Nitrate Deposition 

Two additional effects of a GRB impact, first discussed for general ionizing events by 
Reid et al. (1978), are reduction of sunlight due to opacity in the visible of N0 2 and nitrate 
deposition due to rainout of NO y in the form of HNO3 (nitric acid rain). Here we quantify 
these two effects using our modeling results for the various latitude and time of year cases 
discussed in Sect. 3.2 for a 100 kJ/m 2 fluence burst. A subset of the results in this section 
are discussed elsewhere (Melott et al. 2005). 

3.6.1. Opacity Due to NO 2 

As discussed above, NO2 is a major compound generated by a GRB impact. It has a 
major role in O3 depletion, but also absorbs strongly in the visible, giving it a brown cast. 
Such absorption may easily lower global temperatures, if sufficient NO2 is formed. In order 
to calculate the reduction in transparency due to NO2, we first find the surface irradiance 
at any given time in the presence of NO2. This was done by convolving the solar spectral 
irradiance with the transparency of the Standard US Atmosphere (ASTM Subcommittee 
G3.09. 1999) and the transparency of the N0 2 . The transparency of N0 2 is calculated using 
the Beer-Lambert type relation, similar to that used in this work to calculate UV irradiance, 
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as discussed in Sect. 2.4. The optical depth was computed using the absorption cross section 
of NO2 (Vandaele et al. 1998) and the column density through the atmosphere (computed 
in the model) at the angle between the Sun direction and a normal to the Earth's surface. 
The surface irradiance is then integrated over wavelength and multiplied by the cosine of 
the angle between the Sun and the surface normal to account for insolation. The irradiance 
was averaged over the course of a day. This calculation was repeated with the NO2 column 
density set to zero to obtain the average daily irradiance in the absence of N0 2 . The relative 
transparency was taken as the ratio of these two mean irradiance values, and is shown in 
Figure 19, for the burst cases discussed in Sect. 3.2. We show a full year pre-burst, so that 
the change due to the burst is more obvious. 

The greatest intensity and duration of reduction in sunlight reaching the surface occurs, 
not surprisingly, for those burst cases where the most NO y column density is produced in a 
local area (see Fig. 4). Our results are suggestive that climate change may be possible from 
the several years' reduced sunlight. For example, an estimated reduction of only 0.36% in 
solar flux during the Maunder Minimum may have caused the "Little Ice Age" in Europe 
(Hoyt & Schatten 1993). In our simulations, the greatest NO2 buildup takes place at high 
latitudes, partially persisting during polar summer, which ought to have a greater effect, 
reducing the melting of ice there. A glaciation accompanied the late Ordovician extinction. 
Its initiation is not well understood, but may have required a perturbation such as that 
found here in order to trigger the instability (Herrman et al. 2003, 2004). Further progress 
in understanding this process will require detailed climate modeling. 



3.6.2. Nitrate Deposition 

NO2 reacts with hydroxyl (OH, a product of dissociating water vapor) to make nitric 
acid, HNO3. This is precipitated out in rain or snow, which is one of the primary ways the 
atmosphere returns to equilibrium after our model burst. The global mean surface density 
of nitrate is interesting, of order the amount deposited by lightning (6 x 10~ 3 g/m 2 ) over 
several years, or the amount used in a typical agricultural application. Nitrogen is essential 
for life, but atmospheric N2 is nearly unavailable to most organisms due to the strength 
of the nitrogen triple bond. Most biota respond strongly with increased growth rates due 
to nitrate deposition (Schlesinger 1997). Acids would stress portions of the biosphere (e.g. 
Hatch & Blaustein (2000), but after immediate titration should act as fertilizer (Strong et 
al. 1996; Shi et al. 2005). 

In order to compute nitrate deposition, we used HNO3 rainout data as a reasonable 
measure, since this is the primary pathway by which a GRB would increase nitrate deposition. 
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Rainout is directly computed by the atmospheric model and is empirically based. Using both 
HNO3 concentration and rainout coefficients from the model, we computed HNO3 rainout 
rates as a function of time and latitude. Significant rainout builds up at least a few months 
after the burst and the peak yearly deposition occurs in the third year after the burst in all 
cases presented here. 

Our rainout estimates, and in particular their geographic pattern, must be regarded as 
only exemplary, since they are based on current empirical data. Rainout is more strongly 
affected than are ozone depletion and radiation transparency by the configuration of surface 
features, such as land versus sea, mountains, etc., and overall temperature, which would be 
different at other times. These results should be viewed as a rough guide to what may be 
expected. 

Figure 20 shows nitrate rainout rates (flux of N as nitrate) pre- and post-burst for the 
cases discussed in Sect. 3.2 for a 100 kJ/m 2 fluence burst. Much of the latitude asymmetry 
in this case is due to the fact that there is presently more precipitation in the northern 
hemisphere, especially north of about 45° latitude, as compared to that south of 45° latitude. 
This is primarily due to warmer temperatures at high northern latitudes as compared with 
high southern latitudes (warmer air can hold more moisture) (Stull 2000). More land mass 
at high northern latitudes also contributes, since one trigger of precipitation is uplift flow 
over mountains. 

This latitude asymmetry would not be a feature expected at other times. However, we 
believe our global averages (see Table 5) are reasonable guidelines to effects at earlier times. 
Our baseline rainout rate is empirical and will also include a small additional component 
due to biological sources. N fixation by land plants would not be a significant contribution 
at the Ordovician. However, most biogenic nitrogen is directly deposited on land or in water 
and not from rainout. The background rate due to lightning constitutes at most 70% of 
the pre-burst background as shown in this Figure (Schlesinger 1997). Thus, a burst at the 
Ordovician or earlier would cause a somewhat greater relative change in nitrate rainout rate. 

Table 5 lists global average annual nitrate rainout above baseline for all cases presented 
in Fig. 20. For reference, the baseline (background, without a burst) global average annual 
rainout value is 8.33 x 10~ 3 g/m 2 . The values for several burst cases (especially those 
at +90° and +45°) are similar or somewhat greater than the total generated by lightning 
(6 x 10 -3 g/m 2 ) and other non-biogenic sources (Schlesinger 1997) and values for all cases are 
at least the same order of magnitude. Greater values for bursts at high northern latitudes are 
partly due to the higher levels of rainout in the northern hemisphere in our model. However, 
we also note that for a burst in September over the equator the value is close to that for 
lightning. 
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4. Uncertainties 

Several sources of uncertainty exist in our modeling results. Below we discuss issues as- 
sociated with our computation of ionization profiles; possible effects of redistributed gamma- 
ray photons; possible effects of the GRB afterglow; self-limitation of the production of NO y ; 
the lack of some reactions in the model; the lack of temperature and dynamics feedback in 
the model. 



4-0.3. Ionization Profiles 

Our ionization profiles are computed (as described in Sect. 2.2) using simple energy- 
dependent attenuation coefficients, rather than by a full radiative transfer calculation. We 
have adopted this technique from Gehrels et al. (2003), modifying the functional form of the 
spectrum. The altitude of maximum energy deposition in Gehrels et al. (2003) is around 
34 km, which compares favorably with an altitude of around 32 km from a full radiative 
transfer method (as reported in Smith et al. (2004)). In our case, as shown in Fig. 21, we 
find the altitude of maximum energy deposition to be similar, around 33 km. 

4.0.4. NO y Production 

Some uncertainties in the model's computation of the production of NO y exist, primarily 
due to the lack of some reactions involving NO in the model. The reaction N + N — > N 2 + 7 
is not included in the model and could limit production of NO by reducing the amount of N 
available to react with 2 . However, this reaction is unlikely due to its limited phase space 
and the simultaneous requirement of momentum and energy conservation. A somewhat more 
likely reaction which is not included in the model isN + N + M^N2 + M-|-7, where M is 
an arbitrary atmospheric constituent, effectively either N2 or O2. This is not likely to be a 
large effect, since for our 100 kJ/m 2 fluence burst the N abundance compared to 2 or N 2 
is small, with N/0 2 ~ N/N 2 ~ 10~ 6 . Therefore, it appears much more likely that N would 
react with O2 rather than N in the reaction above. 

There is a limiting reaction for NO production involving excited state N( 2 D) which 
is not included in the model: N( 2 D) + NO -> N 2 + O. The destruction rate of NO with 
excited state N( 2 D) is about twice that of the destruction rate with ground state N( X S). To 
better quantify this effect, we have performed simplified off-line computations to determine 
the maximum production of NO for various ratios of N( 2 D)/N( 4 S). For a 100 kJ/m 2 fluence 
burst, these computations (at 250 K) indicate production concentrations between 1.94 x 10 11 
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and 2.88 x 10 11 cm -3 for ratios between 0.0 and 1.0. This compares favorably with a maximum 
production in the model of 2.5 x 10 11 cm" 3 . This comparison, and the small variation of 
the production concentration with varying N( 2 D)/N( 4 S) ratio indicates that the exclusion of 
this particular reaction is unlikely to be a significant uncertainty. 

Another possible limiting reaction for NO production which is not included in the model 
is H + NO + M — * HNO + M. The rate for this reaction is about 10 -13 cm 3 /s compared to 
10~ n cm 3 /s for the reaction N + NO — > N 2 + O. This could be significant if there is a very 
large concentration of H compared to N. We do note a brief spike in the model results 
where H concentration exceeds N by 2-3 orders of magnitude for a few minutes and then 
quickly drops to 2-3 orders of magnitude smaller than N. Hence, while there is likely a small 
perturbation here which we do not include, it is unlikely to be significant in the long run. 

4-0.5. Additional Radiation Effects 

Our current modeling does not include the redistributed UV from the prompt GRB 
emission. This radiation would be present in the stratosphere where it may affect ozone in 
two competing ways. At less than 242 nm, UV will generally increase ozone by dissociating 
2 . The two oxygen atoms then form ozone by reacting with 2 . UV at greater than 242 
nm will generally decrease ozone because it breaks ozone apart and the oxygen atom has a 
chance of reacting with some other chemical rather than reacting with 2 to reform ozone. 
Full treatment of this issue would require radiative transfer computations which are beyond 
the scope of the present work. However, we have performed a simple test of the magnitude 
of such effects by increasing the solar flux used in the model by 10 4 for 10 s to simulate the 
redistributed UV. The value 10 4 was chosen based on conversations with J. Scalo who has, 
along with D. Smith and J.C. Wheeler, performed radiative transfer computations for cases 
similar to ours Smith et al. (2004). The resulting long-term maximum ozone depletion is 
a few percent greater (42% globally averaged), compared to a case without this increased 
solar flux (36% globally averaged). This test is not definitive since we have not attempted 
to accurately model the spectrum of the redistributed UV. The results at least indicate that 
our results may be somewhat conservative. 

We have not included in our modeling radiation from the afterglow of the GRB. While 
not strictly an uncertainty in the modeling itself, inclusion of afterglow radiation could affect 
our results. The ionizing fluence in afterglow X-ray may be comparable to that of the prompt 
gamma-rays, or somewhat less (Vanderspek et al. 2003; Marshall & Swank 2003). While the 
energy of the photons is lower (0.5 — 2keV), this radiation would arrive over a longer time 
period (days). This would likely lead to an modest increase in the NO y production and 
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O3 depletion, making our current atmospheric results somewhat conservative at worst. A 
perhaps more significant effect is that much of the radiation would arrive after the ozone 
layer has been depleted, at least at some latitudes. This would lead to an enhanced UVB 
flux at the ground, from redistributed X-ray as well as direct UV from the GRB afterglow, 
which would arrive over days to weeks. This would augment the solar UVB, increasing our 
DNA damage values. 

We have also neglected any variability of the parameter values used in the Band spec- 
trum (see Sect. 2.2) or temporal variability of the burst. This would have some effect on 
our ionization profiles which are input to the atmospheric model. We cannot say for certain 
how large of an effect this would be, but we expect it to be small, partly due to the fact that 
our results are not significantly affected by spreading the total energy input over durations 
ranging from 1 s to 1000 s. So far, the total energy appears to be the most important factor 
for long term results. 

Some bursts have values of (3 > — 2 in the Band spectrum, in which case E is not a 
peak of the energy distribution. In these cases more photons are at higher energies. We did 
not include this consideration in our modeling. Long term ozone depletion could be slightly 
enhanced in these cases due to more high energy photons (and so more ionization), but if 
the total energy input is kept the same then we do not expect a large change to our results. 

4-0.6. Temperature and Dynamics Feedback 

The GSFC atmospheric model does not include feedback on temperatures (and thereby, 
dynamics) due to changes in constituent values. This will introduce some uncertainty to 
our results, since depleting ozone changes the energy absorption properties of portions of 
the atmosphere (particularly the stratosphere). This will have complicated effects. Loss of 
ozone leads to cooling which slows depletion reactions involving NO y . However, in the polar 
regions cooling will enhance formation of polar stratospheric clouds which increase ozone 
depletion (see Sect. 3.2). Rosenfield et al. (2002) have investigated the impact on ozone 
depletion by COvmediated stratospheric cooling. For temperature changes of about 5 K 
they find a few percent difference in ozone levels with and without the cooling effect. They 
also find that effects are highly dependent upon latitude and time of year. 

While it is difficult to know how our results may be affected, some estimates may be 
made. Using a ID temperature-constituent coupled model we have determined that changes 
in stratospheric temperature of order 10 K are probable. Key ozone depletion reactions may 
see rate reductions of 30 - 50%. This does not simply translate into changes in ozone depletion 
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however, since there is feedback between ozone levels and temperatures, complicating the 
situation. To get some sense of the effect on ozone levels due to a temperature reduction 
of 10 K, we have run an extreme case where all temperatures everywhere in the model are 
reduced by this amount. This run (in March at the equator, for a 100 kJ/m 2 fluence burst) 
results in maximum localized ozone depletions which are a few (2-3) percent less than for 
a comparable run without the temperature change. Therfore, we conclude that the impact 
of such a temperature change on our primary results (ozone depletion) is not likely to be 
significant. 

Another effect of lack of temperature feedback is on dynamics (transport) in the model. 
Temperatures (and gradients in temperature) affect pressures (and gradients in pressure) 
which drive dynamics. Simple calculations indicate that pressure changes due to a 10 K 
temperature change are relatively minor; at least two orders of magnitude smaller than those 
associated with major weather systems such as tornadoes. However, due to the complications 
of gradients and heating rates, it is not possible to estimate the full effect accurately. The 
primary impact on our results of such uncertainties is likely to be changes in how rapidly 
NO y compounds are dispersed through the atmosphere. This could affect the latitude and 
time distribution of ozone depletion, even if total depletion levels remain similar. 



4.0.7. UVB Calculations 

We have not used a radiative transfer calculation to determine the solar UVB flux at 
the Earth's surface, for use in computing DNA damage. As discussed in Sect. 2.4, we use the 
Beer-Lambert relation and only consider effects of ozone absorption, neglecting scattering 
effects. We believe this is justified for several reasons. First, according to Madronich (1993), 
in the UV band atmospheric absorption (primarily by O3) is dominant over scattering due 
primarily to the relative cross sections for these processes. (The attenuation of a beam by O3 
absorption is approximately two orders of magnitude greater than by scattering, neglecting 
diffuse background.) In addition, while absorption leads to loss of energy, scattering primarily 
changes direction, not energy, leading to a diffuse background of similar flux as would occur 
without scattering. 

We have checked our UVB irradiance results for the unperturbed, equilibrium atmo- 
sphere in comparison with measured values. The typical UVB irradiance at noon in summer 
at the equator today is around 2 — 4 W/m 2 . Using the 1976 U.S. Standard Atmosphere 
ozone column density of 8.1 x 10 18 cm -2 in our code at noon at the summer solstice at the 
equator yields a result of about 3 W/m 2 , in good agreement with observations. Therefore, 
we feel confident that neglecting scattering here does not introduce significant errors. 
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4-0.8. Application to Other Geologic Eras 

As discussed in Appendix A, the GSFC atmospheric model is tuned to the present 
day atmosphere, especially in regard to temperatures and dynamics. Certain aspects of our 
results, such as greater ozone depletion in the southern polar regions and higher rainout 
values in the northern hemisphere, are largely dependent upon continental arrangement. 
Thus, features such as a strong south polar vortex and enhanced PSC formation in that region 
may not be accurate for other periods in the Earth's history, when continental arrangements 
were not the same as today. The latitude dependence of our nitrate deposition results are 
particularly suspect in application to different time periods. 

We do, however, note that latitude dependence of ozone depletion is much more strongly 
affected by differences in what time of year a burst occurs at, as well as over what latitude 
(see Fig. 6). In particular, the North-South asymmetry in ozone depletion which appears for 
a March, equatorial burst is almost completely reversed for a burst over the same latitude 
in September. Since the time of year and latitude of a burst are essentially random, we 
feel reasonably confident that at least our ozone results (as well as DNA damage results, 
which depend upon ozone levels) may be applied to eras when the continental configuration 
differs from today. The random variability of burst latitude and season will overwhelm any 
systematic error due to geography. 

Differences in atmospheric composition over the last billion years could present a prob- 
lem for applying our results. Oxygen levels are the most important factor for our modeling 
since significant 2 must be present in the atmosphere in order to have an ozone shield. The 
oxygen level over the last billion years was generally close to today, except during the Per- 
mian period, and a significant level of ozone should have been in place over this time span 
as well (Berner et al. 2003; Graedel & Crutzen 1993). Carbon dioxide levels have varied 
more greatly and were probably higher at earlier times Kasting (1993); Berner & Kothavala 
(2001), but this is not an important factor for our results. 

Differences in temperatures between today and other eras may introduce more serious 
problems for applying our modeling. Differences in atmospheric temperatures at different 
latitudes affect transport and so may change how quickly and to what extent NO y compounds 
produced by the burst are transported from the latitude over which the burst occurs to other 
locations. This effect would depend especially on temperature differences between the poles 
and tropics, as this affects transport via the "cells" of wind flow in the atmosphere. We 
view our modeling as a "best approximation" in this regard, especially since determining 
atmospheric temperature profiles for past eras is difficult or impossible. 
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5. Discussion 

It has been suggested (Melott et al. 2004) that a GRB may have initiated the second 
largest mass extinction in the Earth's history, that of Late Ordovician period (approximately 
443 My ago). Patterns of extinction at this time are compatible with heightened solar 
UV radiation. Such patterns include dependence on how deep an organism lived in the 
water column and how long an organism spent in a surface-dwelling planktonic larvae state 
(Chatterton & Speyer 1989; Brenchley et al. 1995; Sheehan 2001; Brenchley et al. 2003). In 
addition, surviving fauna which repopulated the globe after the extinction event appear to 
have come from areas of deeper water and higher latitude (Melchin & Mitchell 1991; Sheehan 
2001; Jablonski 2004). Our DNA damage results (Sect. 3.5) indicate that a GRB would most 
likely affect mid to low latitudes and would tend to leave survivors preferentially at higher 
latitudes. 

The Ordovician extinction is associated with a glaciation event, surrounded by periods 
of stable warm climate. Paleoclimate models (Herrman et al. 2003, 2004) have indicated 
that glaciation at this time is not likely to have occurred without some external forcing 
mechanism. Given that even a few lOths of a percent decrease in sunlight reaching the 
Earth's surface can have an effect on global temperatures (Hoyt & Schatten 1993), it seems 
likely that the increase in N0 2 due to a GRB (Sect. 3.6.1) could have provided the necessary 
perturbation to trigger this climate change. 

From our nitrate rainout results (Sect. 3.6.2) we conclude that a burst could supplement 
usual sources of nitrate, and give a modest boost to the transition of plants to land (Gensel 
& Edwards 2001), which was just beginning at the time of the Late Ordovician. Biota 
are extremely responsive to increased nitrate flux (Strong et al. 1996; Schlesinger 1997; Shi 
et al. 2005). Increased productivity is indicated by isotopic excursions in <5 13 C at the late 
Ordovician extinction, coincident with 5 18 indicating lower temperatures (Patzkowsky et 
al. 1997; Sheehan 2001; Brenchley et al. 2003; Pancost et al. 2003). Increased productivity 
could lower atmospheric C0 2 , providing further positive feedback toward glaciation. 

While nitrate deposition could aid land plants, it might at the same time provide further 
stress to some aquatic animals, especially in combination with increased UVB flux. Hatch 
& Blaustein (2000) found significantly reduced survival and activity levels of larval frogs 
in the presence of a combination of increased UVB, nitrate and low pH. This combination 
of stressors could occur, at least temporarily, in shallow water environments in the months 
following a burst, contributing to initiation of the mass extinction. 

We have considered biological effects in the context of mass extinction. However, less 
energetic bursts (or those at greater distances) could also have an interesting impact on 
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mutation rates, and thereby on evolution. Scalo & Wheeler (2002) have considered this 
effect due to the prompt flash of redistributed UV at the surface. Our DNA damage results 
indicate that effects due to increased solar UVB may be just as important. This effect may 
actually be more important, since the duration is much longer (months versus seconds). 

We have not considered the effects of any cosmic rays that may accompany a GRB. 
These would have further effects on atmospheric chemistry (including ozone depletion), would 
produce radioactive nuclides, and also have more direct effects on life through the production 
of showers of penetrating muons (e.g. Thorsett (1995); Dar et al. (1998)). 

We are currently conducting further studies of the effects of including the GRB afterglow 
and varying the spectral parameters used. A related question is the importance of X-ray 
rich GRBs (E ~ 10 keV) and X-ray flashes (E ~ 1 keV) (see, for instance, Zhang et al. 
(2004); Sakamoto et al. (2004); Butler et al. (2005). It is likely that these are GRBs observed 
slightly off-axis. Their galactic population may be larger than that of GRBs, but they are less 
luminous. Hence, these events may present yet another source of astrobiologically interesting 
radiation on a scale similar to GRBs and supernovae. 

As our results indicate, a GRB within 2 kpc of the Earth could cause significant damage 
to the biosphere through ozone depletion, and may have effects on climate through the 
opacity of N0 2 . The "shot" of nitrates that would likely be deposited after a burst could both 
aid land plants in gaining ground and possibly add additional stress to aquatic organisms. 
Such a burst is likely to have occurred within the last billion years and therefore should be 
considered, along with SNe, impacts and other events as contributors to the development of 
life (through both extinctions and enhancements of the mutation rate), both on Earth and 
on other Earth-like planets which may exist elsewhere. 
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A. Description of Atmospheric Model 

For this study we have utilized the NASA Goddard Space Flight Center two dimensional 
photochemical transport model, created and updated by Charles Jackman, Anne Douglass, 
Richard Stolarski, Eric Fleming, and David Considine. The model has been described in 
several publications, including Douglass et al. (1989); Jackman et al. (1990); Considine, 
Douglass, & Jackman (1994); Jackman et al. (1996); Fleming et al. (1999) and Jackman et 
al. (2001). The following description draws upon these articles as well as an unpublished 
overview of the model by Francis Vitt (1994). 

The model contains 65 chemical species (see Table 6). A "family" approach is used for 
the transport of most species. Families include, for instance, O x (0 3 , O, 0( 1 D)), NO y (N, NO, 
N0 2 , N0 3 , N 2 5 , HN0 3 , H0 2 N0 2 , C10N0 2 , BrON0 2 ), HO x (H, OH, H0 2 ), Cl y (CI, CIO, 
HC1, HOC1, C10N0 2 ), etc. It is assumed that family member species are in photochemical 
equilibrium with each other since reactions producing interchanges among members are fast 
compared to photochemical and transport processes that produce net change in the family 
concentration. Several other species are transported separately. See Table 7 for a list of the 
transported families and species. 

Listed in Table 8 are lower (i.e. ground level) boundary conditions for the transported 
species in the model. Species which are solely anthropogenic have been set to zero (both 
boundary conditions and number densities at all altitudes and latitudes). These include 
CFC-11, -12, -22, -113, -114, and -115; Halon-1211, -1301, and -2402; HCFC-141b, -142b, 
and -123. This was done because we are primarily interested in applying our results to 
pre-industrial time periods. 

Some species included have natural as well as anthropogenic sources (such as N 2 0, 
C0 2 , CH 4 , CH3CI, CH 3 Br). We have set these to pre-industrial values (see Table 8). These 
compounds, especially C0 2 , CH 4 and N 2 have probably varied by as much as an order of 
magnitude or more over the last billion years, with C0 2 generally decreasing, CH 4 fluctuating 
and N 2 increasing (see Kasting (1993) and Berner & Kothavala (2001)). The oxygen level 
over the last billion years was generally close to today, except during the Permian period, 
and a significant level of ozone should have been in place over this time span as well (Berner 
et al. 2003; Graedel & Crutzen 1993) 
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A.l. Species Continuity Equation 

The main calculation of the model is the solution of the longitudinally averaged species 
continuity equation: 

Oti ■ 

-^1 = -v ■ Vrii - D{n t ) + S(m) (Al) 

where rij is the species concentration (number density) and v is the wind velocity vector. 
The first term on the right side is advection by circulation. D(rii) is the diffusion term, 
representing small-scale mixing and also effects of planetary and gravity waves. S(rii) is the 
photochemical source term, which is equal to the difference in production and loss terms and 
is dependent upon reaction rate constants, photolysis rates and other species concentrations. 

Under the assumption that changes in species concentration for one time step are small, 
the photochemical terms in the continuity equation can be evaluated using the previous time 
step's values. The species continuity equations are solved by a process splitting method. 
That is, it is assumed that the continuity equation can be separated into a product of 
operators (McRae et al. 1982). It is assumed that advection of species can be separated 
from photochemistry and diffusion. Here, advection operates on the concentration field 
before diffusion and photochemistry. 

The source term, S(rii), is written 

S(m) = P- L(m) (A2) 

where P is the sum of all production sources for species % and L is the loss frequency. For 
species in which the loss process is not linear, L is represented by a Taylor series. The 
continuity equation Al then is written as a finite difference equation: 

n t_ n t-i f l _ n t-i „dL, ,dL. . . „. 

"aT" = + P ~ L(n) + "fln 1 - - " ST 1 - " D(n) (A3) 

where the subscript % has been dropped and n* is the result, n^ 1 is the initial condition 
(from the previous time step), and n is the advected field. 



A. 2. Transport 

The advection scheme includes vertical and meridional (horizontal, north-south) winds. 
These winds are climatological in nature and based on empirical data sets. The primary 
components in deriving the winds are temperature data from the National Centers for En- 
vironmental Prediction (NCEP) and heating rates from climatological distributions of tem- 
perature, ozone, and water vapor. Temperatures and heating rates are updated daily in the 
model and winds are calculated at that time. 
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The advection scheme is described in detail in Fleming et al. (1999). It is mass con- 
serving and uses an upstream piecewise parabolic method of solution (Colella & Woodward 
1984; Carpenter et al. 1990). A time step of 12 hours is used for advection (while the time 
step for constituent changes due to photochemistry is 1 day). The coefficients of the stream 
function equation depend on zonal mean temperature which is based on the 17 year average 
(1979-1995) NCEP temperature data for the troposphere and stratosphere and the CIRA-86 
empirical reference model for the mesosphere above 1 hPa (about 47 km altitude) (Fleming 
et al. 1990). Outside the tropics, winds are derived from temperature using the gradient 
wind relation with measurements from the UARS satellite for tropical zonal winds (20° S - 
20°N). Small scale horizontal mixing (eddy diffusion) coefficients are calculated from me- 
chanical forcing due to gravity waves in the mesosphere and upper stratosphere and from 
the vertical temperature gradient in the troposphere and lower stratosphere. 



A. 3. Photolysis and Chemical Reactions 

The model's binary and tertiary gas-phase-only reactions are listed in Tables 9 and 
10. Photolysis reactions are listed in Table 11. Rates and photolysis cross sections are 
taken from the Jet Propulsion Laboratory (JPL) 2000 recommendations (Sander et al. 2000). 
Reaction rates are updated daily along with temperatures. There are 39 wavelength intervals 
used in the radiative transfer scheme. Multiple scattering is applied and daytime averaged 
photolysis rates are computed every 10 days using a two-stream radiative transfer method 
(Herman 1979). Solar photon fluxes for the wavelength intervals at the top of the atmosphere 
are taken from measurements made by the Solar Stellar Irradiance Comparison Experiment 
(SOLSTICE) and Solar Spectra Irradiance Monitor (SUSIM) instruments on the Upper 
Atmosphere Research Satellite (UARS) (Jackman et al. 1996), and are reproduced here in 
Table 12. 

In the long term version of the model (where daily averages are used), diurnal behavior 
of individual species must be taken into account to calculate photochemical production and 
loss. Species which disappear quickly after sunset, such as 0( 3 P), 0( 1 -D), N, NO, CI, H, 
OH and HO2, are assumed to be zero at night. Daytime average values of NO2 and CIO are 
are needed to calculate odd oxygen loss and these values depend upon the concentrations of 
N2O5 and CIONO2 which build up at night and decrease during the day. Sunset values of 
these are calculated from daytime average values and nighttime production at the expense 
of NO2 is calculated at each grid point for each time step. 

Due to the importance of diurnal behavior, and the short duration of a GRB, we have 
used the short-term version of the atmospheric model (with time step 1 s) to compute the 
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chemistry for several days around the time of the burst. Diurnal behavior is important 
primarily because several constituents important for our study are strongly affected by pho- 
tolysis (see Sect. 3.1 for more discussion of this). 

Important reactions for our purposes in this work occur on timescales of a few to 
tens of seconds. For instance, in the stratosphere NO creation (N + 2 — > NO + O) takes 
about 5 s, NO destruction (N + NO — > N 2 + O) takes about 30 s, and O3 depletion by NO 
(NO + O3 — > NO2 + O2) takes about 85 s. Therefore, our 1 s time step seems adequate to 
accurately model these reactions around the time of the burst. It should be noted that 
destruction of NO proceeds faster with increased concentrations of NO, down to a timescale 
of around a second at the maximum values seen in our results. There is hence a feedback 
between creation and destruction of NO, with more creation leading to faster destruction. 
This will limit the amount of NO created overall and thereby limit ozone depletion as well. 
However, normal diurnal variations of NO (and N0 2 , which also acts to deplete ozone) have 
a strong effect on the total amounts of these compounds which are present at any one time 
(see Sect. 3.1). 



A. 4. Heterogeneous Processes and Galactic Cosmic Rays 

Heterogeneous processes which occur on the stratospheric sulfate aerosol (SSA) layer 
and on polar stratospheric clouds (PSCs) are included as described in Considine, Douglass, & 
Jackman (1994). SSAs are primarily liquid droplets of H2SO4 and H2O. PSCs comprise two 
classes; Type 1 is primarily frozen nitric acid trihydrate (NAT), Type 2 is frozen H 2 0. The 
formation temperatures of PSCs are low and they occur most frequently in the polar regions, 
while SSAs are more uniformly distributed (with latitude) in the stratosphere. Heterogeneous 
processes occurring on the surfaces of these particulates (especially PSCs) plays an important 
role in the depletion of ozone and are partly responsible for the preferential depletion in the 
south polar region where formation of PSCs is more probable. The formation of PSCs 
depends on the concentrations of NO y (particularly HN0 3 ) and H 2 0. 

In the model, computed concentrations of NO y and H 2 are combined with climato- 
logical temperature distributions (obtained from National Meteorological Center data) to 
calculate: the probability of cloud occurrence; the amount of H 2 and HN0 3 removed from 
gas phase in producing clouds; and the surface area density of clouds. Surface area density is 
important in computing rates of reactions which occur on the particulates. It is important to 
note that this is not a microphysical model of PSCs. While the computed characteristics of 
PSCs depends on several assumptions, the results are reasonably insensitive to large changes 
in the assumptions (Considine, Douglass, & Jackman 1994). 
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The effects of galactic cosmic rays (GCRs) are also included in the model. Ionization 
rate profiles produced by GCRs were calculated by Nicolet (1975) for both solar maximum 
and minimum. Rates in the model for a given year are calculated using the yearly averaged 
monthly sunspot number by linear interpolation of Nicolet 's rates. 
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Fig. 1. — Biological weighting function of Setlow (1974) and Smith et al. (1980), quantifying 
damage to DNA by UVB. 
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Fig. 2. — Annual variation of the column density of O3 in an unperturbed run, with scales 
for both Dobson units (left) and 10 18 cm -2 (right). The horizontal axis labels give month 
initials starting in January. (See color figure online.) 
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Fig. 3. — Column density of NO, NO2, NO y , and O3, normalized to their maximum values 
over this time period. These are hourly results from the short-term version of the model for 
a 100 kJ/m 2 burst over the equator in late March at noon. Sudden changes (particularly in 
NO and N0 2 ) occur at sunrise/sunset. (Burst occurs at time 0, which is local noon.) 
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Fig. 4. — Column density of NO y in units of 10 16 cm -2 , for 100 kJ/m 2 bursts over latitudes 
+90°, +45°, the equator, —45°, and —90°, at the equinoxes and solstices. (Bursts occur at 
month 0. See color figure online.) 
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Fig. 5. — Column density of O3 with scales for both Dobson units (left) and 10 18 cm 2 
(right), for 100 kJ/m 2 bursts over latitudes +90°, +45°, the equator, -45°, and -90°, at the 
equinoxes and solstices. (Bursts occur at month 0. See color figure online.) 
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Fig. 6. — Pointwise percent change in column density of O3 (comparing runs with and 
without burst), for 100 kJ/m 2 bursts over latitudes +90°, +45°, the equator, —45°, and 
—90°, at the equinoxes and solstices. (Bursts occur at month 0.) Note that in this plot 
white indicates values equal to or greater than 0.0. The "wedges" of white that appear in 
the polar regions in some frames are the only areas of positive percent change. See color 
figure online. 
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Fig. 7. — Percent change in globally averaged column density of 3 , for 100 kJ/m 2 bursts 
over latitudes +90°, +45°, the equator, —45°, and —90°, at the equinoxes and solstices. 
(Bursts occur at month 0.) 
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Table 1. Maximum globally averaged 3 percent change, for 100 kJ/m 2 bursts. 



Latitude (deg) 


March 


June 


September 


December 


+90 


-26 


-29 


-32 


-31 


+45 


-30 


-32 


-33 


-31 





-36 


-37 


-38 


-36 


-45 


-32 


-31 


-31 


-31 


-90 


-31 


-30 


-26 


-36 
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Table 2. Maximum localized 3 percent change, for 100 kJ/m 2 bursts. 



Latitude (deg) 


March 


June 


September 


December 


+90 


-60 


-66 


-74 


-70 


+45 


-56 


-63 


-72 


-66 





-55 


-55 


-60 


-55 


-45 


-67 


-63 


-55 


-66 


-90 


-71 


-66 


-59 


-62 
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Fig. 8. — Pointwise percent change in number density of NO y (comparing runs with and 
without burst), for 100 kJ/m 2 , March, equator case, 12 months after burst. (See color figure 
online.) 
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Latitude 

Fig. 9. — Pointwise percent change in number density of O3 (comparing runs with and 
without burst), for 100 kJ/m 2 , March, equator case, 12 months after burst. Black in this 
figure represents both positive and negative percent changes. White is set to 0.0. The white 
"boundary" starting at about 15 km in the polar regions and moving upward toward the 
equator separates areas of negative percent change at higher altitudes and positive percent 
change at lower altitudes. See color figure online. 
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Fig. 10. — Intensity of sunlight at noon (in arbitrary units), beginning at the four times 
of year at which we have input bursts. On the upper axes, M = March, J = June, S = 
September, D = December. 
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Table 3. Maximum globally averaged O3 percent change, every 2 hours, for 100 kJ/m 2 , 

March, equator burst case. 



Time of Day 


Global Average 


Localized 


Noon 


-36 


-55 


2pm 


-37 


-55 


4pm 


-37 


-55 


6pm 


-36 


-55 


8pm 


-40 


-57 


10pm 


-40 


-57 


Midnight 


-40 


-57 


2 am 


-40 


-57 


4am 


-40 


-57 


6 am 


-37 


-55 


8 am 


-37 


-55 


10am 


-37 


-55 




Fig. 11. — Column density of NO y in units of 10 16 cm 2 . Hourly results for a 100 kJ/m 
burst over the equator in late March. (Burst occurs at time 0. See color figure online.) 
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Fig. 12. — Column density of O3 with scales for both Dobson units (left) and 10 18 cm -2 
(right). Hourly results for a 100 kJ/m 2 burst over the equator in late March. (Burst occurs 
at time 0. See color figure online.) 
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Months 

Fig. 13. — Column density of NO y in units of 10 16 cm -2 , for 10 kJ/m 2 , March, Equator (top 
panel) and September, +45° (bottom panel) cases. (Bursts occur at month 0.) Note that 
the scale here is not the same as that for Fig. 4. (See color figure online.) 
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Fig. 14. — Column density of O3 with scales for both Dobson units (left) and 10 18 cm -2 
(right), for 10 kJ/m 2 , March, Equator (top panel) and September, +45° (bottom panel) 
cases. (Bursts occur at month 0.) Note that the scale here is not the same as that for Fig. 5. 
(See color figure online.) 
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Fig. 15. — Pointwise percent change in column density of O3 (comparing runs with and 
without burst), for 10 kJ/m 2 , March, Equator (top panel) and September, +45° (bottom 
panel) cases. Note that the scale here is not the same as that for Fig. 6. In contrast to 
Fig. 6, there are no areas of positive percent change in these cases. (See color figure online.) 
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Fig. 16. — Percent change in globally averaged column density of O3, for 10 kJ/m 2 , March, 
Equator (solid line) and September, +45° (dotted line) case. 
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Table 4. Maximum O3 percent changes for 10 kJ/m 2 cases. 



Burst Case 


Global Average 


Localized 


March, Equator 


-16 


-28 


September, +45° 


-13 


-38 
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Fig. 17. — Relative DNA damage (dimensionless), normalized by the annual global average 
damage in the absence of a GRB, for 100 kJ/m 2 bursts over latitudes +90°, +45°, the 
equator, —45°, and —90°, at the equinoxes and solstices. (Bursts occur at month 0.) Note 
that white indicates values between 1.0 and 0.0 (see color figure online). Experiments suggest 
significant mortality for marine microorganisms should certainly exist in areas where this 
measure exceeds 2. 
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Fig. 18. — Relative DNA damage (dimensionless), normalized by the annual global average 
damage in the absence of a GRB, for 10 kJ/m 2 , March, Equator (top panel) and September, 
+45° (bottom panel) cases. (Bursts occur at month 0.) Note that white indicates values 
between 1.0 and 0.0 (see color figure online). Note that the scale here is not the same as 
that for Fig. 17. 
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Fig. 19. — The time development of relative solar fluence reaching the surface for 100 kJ/m 2 
bursts over latitudes +90°, +45°, the equator, —45°, and —90°, at the equinoxes and solstices. 
Relative fluence is the computed amount divided by that in the absence of NO2 absorption. 
The value is set to 1 for total polar darkness. 
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Fig. 20. — Rainout of N as nitrate (NOj) in units of 10~ 10 g/m 2 /s for 100 kJ/m 2 bursts 
over latitudes +90°, +45°, the equator, —45°, and —90°, at the equinoxes and solstices. The 
pre-burst background includes all non-anthropogenic rainout at modern times, and would 
be somewhat lower prior to the advent of terrestrial nitrogen-fixing plants. 
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Table 5. Annual global average nitrate deposition (10 3 g/m 2 ) above baseline in year 3 

(post-burst), for 100 kJ/m 2 bursts. 



Latitude (deg) 


March 


June 


September 


December 


+90 


5.19 


6.52 


10.9 


8.98 


+45 


4.45 


5.62 


9.36 


6.87 





3.52 


3.37 


5.67 


4.17 


-45 


3.83 


2.72 


2.84 


2.72 


-90 


4.40 


2.95 


2.49 


2.68 
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Fig. 21. — Energy deposition (10 5 MeV/cm 3 ) for a 100 kJ/m 2 burst over the equator, in late 
March. 
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Table 6. Species Included in the Model 



Number Species Number Species 



1 


0( 3 P) 


34 


H0 2 N0 2 


2 


OCD) 


35 


Cl y 


3 


o 2 


36 


CFCI3 


4 


o 3 


37 


CF 2 C1 2 


5 


NO 


38 


CCI4 


6 


N0 2 


39 


CH3CI 


7 


N0 3 


40 


liquid H 2 


8 


N 2 5 


41 


o x 


9 


N 


42 


CH3CCI3 


10 


N 2 


43 


2 ( 1 A) 


11 


HN0 3 


44 


N( 2 D) 


12 


N 2 


45 


co 2 


13 


H 


46 


BrO 


14 


OH 


47 


Br 


15 


H0 2 


48 


HBr 


16 


H 2 


49 


BrON0 2 


17 


H 2 2 


50 


Br x 


18 


HC1 


51 


CH 3 Br 


19 


H 2 


52 


CBrF 3 


20 


CH 4 


53 


CBrClF 2 


21 


CO 


54 


CHC1F 2 


22 


CH 3 


55 


C 2 Cl3F 3 


23 


CH 3 2 


56 


C 2 C1 2 F 4 


24 


CH3O 


57 


C 2 C1F 5 


25 


H 2 CO 


58 


HF 


26 


HCO 


59 


CC1FO 


27 


HOC1 


60 


CF 2 


28 


CH3OOH 


61 


BrCl 


29 


CI 


62 


C1 2 2 


30 


CIO 


63 


C10 x 


31 


HC1 


64 


ClONO 
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Table 6 — Continued 



Number 


Species 


Number 


Species 


32 


C10N0 2 


65 


Cl 2 


33 


NO y 
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Table 7. Transported Species and Families 



Number 


Species or Family 


Number 


Species or Family 


1 


o x 


20 


C 2 C1F 5 


2 


NO y 


21 


CBrClF 2 


3 


Cly 


22 


CBrF 3 


4 


N 2 


23 


HF 


5 


CFC1 3 


24 


CC1FO 


6 


CF 2 C1 2 


25 


CF 2 


7 


CCI4 


26 


C10N0 2 


8 


CH3CI 


27 


N 2 5 


9 


CH 4 


28 


N 2 


10 


H 2 


29 


HNO3 (solid) 


11 


CO 


30 


H 2 (solid) 


12 


CH X 


31 


co 2 


13 


CH3CCI3 


32 


CH3CC1 2 F 


14 


HNO3 


33 


CH 3 CC1F 2 


15 


Br y 


34 


CHC1 2 CF3 


16 


CH 3 Br 


35 


C 2 Br 2 F 4 


17 


CHC1F 2 


36 


H 2 2 


18 


C2CI3F3 


37 


HO x 


19 


C2CI2F4 






Note. - 


Families are those with "x" or 


"y" subscript. 
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Table 8. Lower Boundary Conditions for Transported Species 



Species Boundary Condition Type Value 



N 2 


mixing ratio 


270 ppbv 


CH 4 


mixing ratio 


700 ppbv 


co 2 


mixing ratio 


280 ppmv 


CO 


mixing ratio 


100 ppbv 


CH 3 Br 


mixing ratio 


7.0 pptv 


CH3CI 


mixing ratio 


483 pptv 


H 2 


mixing ratio 


500 ppbv 


NO y - HNO3 


mixing ratio 


105 pptv 


CCI4 


mixing ratio 


0.0 


CH3CCI3 


mixing ratio 


0.0 


Cly 


flux 


0.0 cm~ 2 s~ 


Br y 


flux 


0.0 cm~ 2 s~ 


HF 


flux 


0.0 cm~ 2 s~ 


CC1FO 


flux 


0.0 cm~ 2 s~ 


CF 2 


flux 


0.0 cm~ 2 s~ 


CH3OOH 


flux 


0.0 cm~ 2 s~ 


o x 


deposition velocity 


0.1 cm s _1 
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Table 9. Binary Reactions and Rate Constants 



Reaction 



O + 3 -> 2 + 2 
OH + 3 -> H0 2 + 2 
H0 2 + 3 -> OH + 20 2 
CIO + H0 2 -> HOC1 + 2 
CI + H 2 2 -> HC1 + H0 2 
0( 1 D) + M -> 0( 3 P) + M 

NO + 3 -> N0 2 + 2 

N0 2 + 3 -> N0 3 + 2 

H + 3 -> OH + 2 

OH + C10N0 2 -> HOC1 + N0 3 

CH 4 + OH -> CH 3 + H 2 

CH 3 2 + NO -> CH3O + N0 2 

OH + CH3CI -> H 2 + CH 2 C1 

CH3O + 2 -> CH 2 + H0 2 

H0 2 + H0 2 -> H 2 2 + 2 

N + 2 -> NO + O 

CH 2 + O -> HCO + OH 

CH 3 2 + H0 2 -> CH3OOH + 2 

CI + H 2 -> HC1 + H 

CI + 3 -> CIO + 2 

CIO + O -> CI + 2 

CI + CH 4 -> HC1 + CH 3 

HC1 + OH -> CI + H 2 

CIO + NO -> CI + N0 2 

OH + H 2 2 -> H 2 + H0 2 

H 2 + OH -> H 2 + H 

N 2 5 + M -> N0 2 + N0 3 + M 

+ H 2 2 -> OH + H0 2 

O + C10N0 2 -> CIO + N0 3 

CO + OH -> C0 2 + H 



Rate Constant 



ki = 8.0(-12) exp(-2060/T) 
k 2 = 1.5(-12) exp(-880/T) 
k 3 = 2.0(-14) exp(-680/T) 
k 4 = 4.8(-13) exp(700/T) 
k 5 = l.l(-ll) exp(-980/T) 
k 6 = 0.78[1.8(— 11) exp(110/T)] 



+0.21(3.2 


(-11 


) exp(70/T)] 


k 7 = 


= 3.0(- 


-12) 


exp(- 


-1500/T) 


k 8 = 


= l-2(- 


-13) exp(- 


-2450/T) 


kg - 


= l-4(- 


-10) exp(- 


-470/T) 


kio 


= 1.2( 


-12) 


expi 


-330/T) 


kn 


= 2.45(-12) exp(-1775/T) 


ki 2 


= 3.0( 


-12) 


expi 


280/T) 


ki 3 


= 4.0( 


-12) 


expi 


-1400/T) 


ki4 


= 3.9( 


-14) 


expi 


-900/T) 


ki 5 


= 2.3( 


-13) 


expi 


600/T) 


ki 6 


= 1.5( 


-11) 


expi 


-3600/T) 


ki 7 


= 3.4( 


-11) 


expi 


-1600/T) 


ki 8 


= 3.8( 


-13) 


expi 


800/T) 


ki 9 


= 3.7( 


-11) 


expi 


-2300/T) 


k 20 


= 2.3( 


-11) 


expi 


-200/T) 


k 2 i 


= 3.0( 


-11) 


expi 


-70/T) 


k 22 


= 9.6( 


-12) 


expi 


-1360/T) 


k 23 


= 2.6( 


-12) 


expi 


-350/T) 


k 24 


= 6.4( 


-12) 


expi 


290/T) 


k 25 


= 2.9( 


-12) 


expi 


-160/T) 


k 26 


= 5.5( 


-12) 


expi 


-2000/T) 


k 27 


= 2.7( 


-27) 


expi 


-11000/T) 


k 28 


= 1.4( 


-12) 


expi 


-2000/T) 


k 29 


= 2.9( 


-12) 


expi 


-800/T) 


k 3 o 


= 1.5( 


-13) 







Table 9 — Continued 



Reaction 



HN0 3 + OH -> N0 3 + H 2 



NO + H0 2 -> OH + N0 2 

H 2 + 0( 1 D) ^OH + OH 

OH + H0 2 -> H 2 + 2 

OH + O -> H + 2 

H0 2 + O -> OH + 2 

N0 2 + O -> NO + 2 

N 2 + 0( 1 D) ^NO + NO 

N + NO -> N 2 + O 

H 2 + 0( 1 D) ^OH + H 

CH 4 + 0( 1 D) -> CH 3 + OH 

H 2 CO + OH -> H 2 + HCO 

HCO + 2 -> CO + H0 2 

CI + H0 2 -> HC1 + 2 

CCI4 + 0( 1 D) -> 4C1 + Prods 

OH + H0 2 N0 2 -> H 2 + 2 + N0 2 

CH 4 + 0( 1 D) -> H 2 + H 2 CO 

OH + CH3OOH -> H 2 + CH 3 2 

OH + OH^H 2 + 

CIO + OH -> CI + H0 2 

CIO + OH -> HC1 + 2 

HOC1 + OH -> H 2 + CIO 

Cl + H 2 CO ^HCl + HCO 

H0 2 + H0 2 + M -> H 2 2 + 2 + M 

CFC10 + 0( 1 D) -> Prods 

CF 2 + 0( 1 D) -> Prods 

CI + H0 2 -> OH + CIO 

N + OH^NO + H 



Rate Constant 



k 31 = 7.2(-15) exp(-785/T) 
+(N d 1.9(-33) exp(-725/T)/ 
(1 + N d 1.9(-33) exp(-725/T)/ 



(4.1(- 


-16) 


exp( 


-1440/T)))) 


k 32 = 


3.5( 


-12) 


exp(250/T) 


k 3 3 = 


2.2( 


-10) 




k 34 = 


4.8( 


-11) 


exp(250/T) 


k 35 = 


2.2( 


-11) 


exp(120/T) 


k 3 e = 


3.0( 


-11) 


exp(200/T) 


k 3 7 = 


5.6( 


-12) 


exp(180/T) 


k 38 = 


6.7( 


-11) 




k 39 = 


2.1( 


-11) 


exp(100/T) 


k 40 = 


1.1( 


-10) 




k 41 = 


1.15 


>(-10) 


k 42 = 


1.0( 


-11) 




k 43 = 


3.5( 


-12) 


exp(140/T) 


k 44 = 


1.8( 


-11) 


exp(170/T) 


k 45 = 


3.3( 


-10) 




k 46 = 


1.3( 


-12) 


exp(380/T) 


k 47 = 


8.0( 


-12) 




k 48 = 


3.8( 


-12) 


exp(200/T) 


k 49 = 


4.2( 


-12) 


exp(-240/T) 


k 5 o = 


7.4( 


-12) 


exp(270/T) 


k 5 i = 


3.2( 


-13) 


exp(320/T) 


k 5 2 = 


3.0( 


-12) 


exp(-500/T) 


k 53 = 


8.1( 


-11) 


exp(-30/T) 


k 54 = 


1.7( 


-33) 


exp(1000/T) 


k 5 5 = 


1.9( 


-10) 




k 5 e = 


7.4( 


-11) 




k 5 7 = 


4.1( 


-11) 


exp(-450/T) 


k 5 s = 


5.0( 


-11) 
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Table 9 — Continued 



Reaction 



Rate Constant 



BrO + NO — > N0 2 + Br 

H0 2 N0 2 + M -> H0 2 + N0 2 + M 

H + H0 2 -> H 2 + 2 

H + H0 2 ^H 2 + 

H + H0 2 -> OH + OH 

NO + N0 3 -> N0 2 + N0 2 

OH + CH3CCI3 -> 3(C1) 

N( 2 D) + 2 -> NO + O 

N 2 + 0( 1 D) -> N 2 + 2 

0( 1 D) + CF 2 C1 2 -> CIO + CI + Prods 

N + N0 2 -> N 2 + O 

0( 1 D) + CFClg -> CIO + 2C1 + Prods 

O^D) + DECAY -> O + hz/ 

2 ( 1 A) + DECAY -> 2 + hi/ 

N( 2 D) + -> N + 

2 ( 1 A) + 2 ^0 2 + 2 

2 ( 1 A) + 3 ^ 2 + 2 + 

2 ( 1 A) + N0 -> 2 + N0 

+ + M^0 2 + M 

Br + 3 -> BrO + 2 

Br + H0 2 -> HBr + 2 

BrO + CIO -> Br + C100 

BrO + BrO -> Br + Br + 2 

OH + HBr -> H 2 + Br 

CH 3 Br + OH -> Br + Prods 

CHC1F 2 + OH -> CI + 2F + Prods 

C 2 C1 3 F 3 + 0( 1 D) -> 3C1 + 3F + Prods 

C 2 C1 2 F 4 + O^D) -> 2C1 + 4F + Prods 

C 2 C1F 5 + 0( 1 D) -> CI + 5F + Prods 

BrO + CIO -> Br + 0C10 

BrO + CIO -> BrCl + 2 



k 5 g 


— 


8.8(- 


-12) exp(260/T) 


k 6 o 


— 


2.1(- 


-27) exp(-10900/T) 


k 6 i 


— 


1.06( 


-11) 


k 62 


— 


1.60( 


-12) 


k 6 3 


= 


6.88( 


-11) 


k 6 4 


= 


1.5(- 


-11) exp(170/T) 


k 6 5 




1.8(- 


■12) exp(-1550/T) 


k 6 e 




5.0(- 


■11) 


k 6 7 




4.9(- 


11) 


k 6 s 




1.4(- 


■10) 


k 6 g 




5.8(- 


-12) exp(220/T) 


k 7 o 




2.3(- 


■10) 


k 71 




6.8(- 


-09) 


k 7 2 


= 


2.6(- 


-10) 


k 7 3 


= 


1.4(- 


10) 


k 74 


— 


3.6(- 


-18) exp(-220/T) 


k 7 5 




5.2(- 


-11) exp(-2840/T) 


k 7 e 




4.5(- 


17) 


k 77 




1.4(- 


•33) exp(408/T) 


k 78 




1.7(- 


■11) exp(-800/T) 


k 79 




1.5(- 


■11) exp(-600/T) 


k 8 o 




2.3(- 


■12) exp(260/T) 


k 8 i 




1.5(- 


-12) exp(230/T) 


k 82 




1.1(- 


-11) 


k 83 




4.0(- 


■12) exp(-1470/T) 


k 84 




1.0(- 


-12) exp(-1600/T) 


k 85 




2.00( 


-10) 


k 86 




1.00( 


-10) 


k 87 




5.00( 


-11) 


k 88 




9.5(- 


-13) exp(550/T) 


k 89 




4.1(- 


■13) exp(290/T) 
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Table 9 — Continued 



Reaction 



C1 2 2 + M -> CIO + CIO 
BrO + O -> Br + 2 
BrO + H0 2 -> HOBr + 2 
Br + CH 2 -> HBr + CHO 
CH 4 + 0( 1 D) -> H + CH 3 
0( 1 D) + CClBrF 2 -> Prods 
0( 1 D) + CBrF 3 -> Prods 
O^D) + CH 3 CC1 2 F -> Prods 
OH + CH3CCI2F -> Prods + H 2 
CI + CH 3 CC1 2 F -> Prods + HC1 
0( 1 D) + CH 3 CF 2 C1 -> Prods 
OH + CH3CF2CI -> Prods + H 2 
CI + CH 3 CF 2 C1 -> Prods + HC1 
0( 1 D) + CHC1 2 CF 3 -> Prods 
OH + CHC1 2 CF 3 -> Prods + H 2 
CI + CHC1 2 CF 3 -> Prods + HC1 
0( 1 D) + C 2 Br 2 F 4 -> Prods 
O^D) + CH 3 Br -> Br + CH 3 
O^D) + CHC1F 2 -> CI + CF 2 + H 
O + CH 3 -> CH 2 + H 
3 + CH 3 -> CH3O + 2 
CH 3 2 + CH 3 2 -> 2CH 3 + 2 
CFC1 3 + OH -> 3C1 + HF + CO 
CF 2 C1 2 + OH -> CI + HF + CC1FO 
CH 3 C1 + CI -> 2C1 + CH 3 



Rate Constant 



kgo = 


= 1.3(- 


-27) exp(-8744/T) 


kgi = 


= 1.9(- 


-11) exp(230/T) 


kg 2 = 


= 3.4(- 


-12) exp(540/T) 


kg3 = 


= l-7(- 


-11) 


cxp(-800/T) 


k g4 = 


= 3.0(- 


-11) 




kg5 = 


= 1.5(- 


-10) 




kg6 = 


= 1.0(- 


-10) 




kg7 = 


= 2.6(- 


-10) 




kgs = 


= l-7(- 


-12) exp(-1700/T) 


kgg - 


= 1.8(- 


-12) exp(-2000/T) 


kioo 


= 2.2( 


-10) 




kioi 


= 1.3( 


-12) 


exp(-1800/T) 


kl02 


= 1.4( 


-12) 


exp(-2420/T) 


kl03 


= 2.0( 


-10) 




kl04 


= 7.0( 


-13) 


exp(-900/T) 


kl05 


= 4.4( 


-12) 


exp(-1750/T) 


kl06 


= 1.6( 


-10) 




kl07 


= 1.8( 


-10) 




kios 


= 1.0( 


-10) 




kiog 


= 1.1( 


-10) 




kno 


= 5.4( 


-12) 


exp(-220/T) 


km 


= 2.5( 


-13) 


exp(190/T) 


km 


= 1.0( 


-12) 


exp(-3700/T) 


kii 3 


= 1.0( 


-12) 


exp(-3600/T) 


kiu 


= 3.2( 


-11) 


exp(-1250/T) 



Note. — Read "Prods" as "Products." Read 1.0(-10) as 1.0 x lO" 10 . 
M = N 2 , 2 (arbitrary third body). N d = total number density at a 
particular grid point. 
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Table 10. Tertiary Reactions and Rate Constants 



Reaction 



£,300 
tv 



n 



1,300 



m 



+ 2 + M^0 3 + M 
H + 2 + M -> H0 2 + M 
OH + OH + M -> H 2 2 + M 
OH + N0 2 + M -> HN0 3 + M 
CIO + N0 2 + M -> C10N0 2 + M 
H0 2 + N0 2 + M -> H0 2 N0 2 + M 
N0 2 + O + M -> N0 3 + M 
N0 2 + N0 3 + M -> N 2 5 + M 
CH 3 + 2 + M -> CH 3 2 + M 
NO + O + M -> N0 2 + M 
0( 1 D) + N 2 + M -> N 2 + M 
BrO + N0 2 + M -> BrON0 2 + M 
CIO + CIO + M -> C1 2 2 + M 
CO + O + M -> C0 2 + M 



6.0( 


-34) 


2.3 


0.0 


0.0 


5.7( 


-32) 


1.6 


7.5(-ll) 


0.0 


6.2( 


-31) 


1.0 


2.6(— 11) 


0.0 


2.5( 


-30) 


4.4 


1.6(— 11) 


1.7 


1.8( 


-31) 


3.4 


1.5(— 11) 


1.9 


1.8( 


-31) 


3.2 


4.7(-12) 


1.4 


9.0( 


-32) 


2.0 


2.2(-ll) 


0.0 


2.2( 


-30) 


3.9 


1.5(— 12) 


0.7 


4.5( 


-31) 


3.0 


1.8(— 12) 


1.7 


9.0( 


-32) 


1.5 


3.0(-ll) 


0.0 


3.5( 


-37) 


0.6 


0.0 


0.0 


5.2( 


-31) 


3.2 


6.9(-12) 


2.9 


2.2( 


-32) 


3.1 


3.5(-12) 


1.0 


2.0( 


-37) 


0.0 


0.0 


0.0 



Note. — Read 1.0(-10) as 1.0 x 10~ 10 . M = N 2 , 2 (arbitrary third 
body) 

Rate constants are given by: 

k = [k (T)]/[l + fc (T)[M]/A; oo (T)]0.6{ 1+ [' Ofll0 ( fe0 ^[ M ]/^( T »] 2 >" 1 

Low pressure limit: k (T) = k™°(T/300)- n (cm 6 s" 1 ) 
High pressure limit: fc 00 (T) = £;^ (T/300)- m (cm 3 s" 1 ) 



Table 11. Photodissociations 



Reaction 



2 -> O + O 
3 ^0 2 ( 1 A) + 0( 1 D) 

3 -> o 2 + o 

H 2 -> H + OH 

N0 3 -> N0 2 + O 

HNO3 -> OH + N0 2 

N0 2 -> NO + O 

H 2 2 -> OH + OH 

N 2 5 -> N0 2 + N0 3 

H 2 CO -> HCO + H 

H 2 CO -> H 2 + CO 

C0 2 -> CO + O 

CH3OOH -> CH3O + OH 

N 2 -> N 2 + O 

C10N0 2 -> CI + N0 3 

NO -> N + O 

N0 3 -> NO + 2 

HC1 -> H + CI 

CCI4 -> 4C1 + Fragment 

CH3CI -> CH 3 + CI 

CFCI3 -> 3C1 + Fragment 

CF 2 C1 2 -> 2C1 + Fragment 

HOC1 -> OH + CI 

H0 2 N0 2 -> OH + N0 3 

H 2 -> H 2 + 0( 1 D) 

CH3CCI3 -> 3C1 + Fragment 

BrO -> Br + O 

BrON0 2 -> Br + N0 3 

CH 3 Br -> CH 3 + Br 

CF 3 Br -> Br + 3F + Fragment 

CF 2 ClBr -> Br + CI + 2F + Fragment 



Table 11 — Continued 



Reaction 



CHC1F 2 -> CI + 2F + Fragment 
C2CI3F3 -> 3C1 + 3F + Fragment 
C 2 C1 2 F 4 -> 2C1 + 4F + Fragment 
C 2 C1F 5 -> CI + 5F + Fragment 
C1 2 2 -> CI + OCIO 
OCIO -> CI + 2 
Cl 2 -> CI + CI 
ClONO -> CIO + NO 
BrCl -> Br + CI 
C0 2 -> CO + O^D) 
H0 2 N0 2 -> H0 2 + N0 2 
CC1FO -> CI + F + Fragment 
CF 2 -> 2F + Fragment 
N 2 -> NO + N 
2 -> O + O^D) 
H0 2 -> O + OH 

CH 3 CFC1 2 — > 2C1 + F + Fragment 
CH 3 CF 2 C1 -> CI + 2F + Fragment 
CF 3 CHC1 2 -> 2C1 + 3F + Fragment 
CF3CHFCI -> 4F + CI + Fragment 
CF 3 CF 2 CHC1 2 -> 5F + 2C1 + Fragment 
CF 2 C1CF 2 CHFC1 -> 5F + 2C1 + Fragment 
CF 2 Br 2 -> 2Br + 2F + Fragment 
CF 2 BrCF 2 Br -> 3Br + 3F + Fragment 
CHBr 3 — > 3Br + Fragment 
CC1 2 -> 2C1 + Fragment 
C10N0 2 -> CIO + N0 2 
CH 4 — > CH 3 + H 
CH4 — > CH 2 + H 2 
CH 4 — > CH + H + H 2 
HOBr -> HO + Br 
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Table 11 — Continued 



Reaction 



CH 3 2 -> CH 3 + 2 

BrON0 2 -> BrO + N0 2 

CH3CFCI2 -> C0 2 + CC1FO + CI + H + H 2 

CH 3 CF 2 C1 -> C0 2 + CF 2 + CI + H + H 2 

CF 3 CHC1 2 -> CC1FO + CF 2 + CI + H 

C 2 F 4 Br 2 -> 2CF 2 + 2Br 
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Table 12. Wavelength Bands and Corresponding Solar Fluxes 



Number Band (nm) Flux (photons cm 2 s r ) 



1 


121.067- 


122.067 


4.006 x 


10 11 


2 


170.0 - 


172.4 


1.764 x 


10 11 


3 


172.4 - 


173.9 


1.017 x 


10 11 


4 


173.9 - 


175.4 


1.302 x 


10 11 


5 


175.4 - 


177.0 


1.722 x 


10 11 


6 


177.0 - 


178.6 


2.200 x 


10 11 


7 


178.6 - 


180.2 


2.438 x 


10 11 


8 


180.2 - 


181.8 


3.215 x 


10 11 


9 


181.8 - 


183.5 


3.691 x 


10 11 


10 


183.5 - 


185.2 


3.549 x 


10 11 


11 


185.2 - 


186.9 


4.190 x 


10 11 


12 


186.9 - 


188.7 


5.555 x 


10 11 


13 


188.7- 


190.5 


6.464 x 


10 11 


14 


190.5 - 


192.3 


7.342 x 


10 11 


15 


192.3 - 


194.2 


7.709 x 


10 11 


16 


194.2 - 


196.1 


1.055 x 


10 12 


17 


196.1 - 


198.0 


1.189 x 


10 12 


18 


198.0 - 


200.0 


1.333 x 


10 12 


19 


200.0 - 


202.0 


1.610 x 


10 12 


20 


202.0 - 


210.5 


1.195 x 


10 13 


21 


210.5 - 


219.8 


3.720 x 


10 13 


22 


219.8 - 


229.9 


5.643 x 


10 13 


23 


229.9 - 


241.0 


6.210 x 


10 13 


24 


241.0 - 


253.2 


8.321 x 


10 13 


25 


253.2 - 


266.7 


2.363 x 


10 14 


26 


266.7- 


281.7 


4.169 x 


10 14 


27 


281.7- 


285.7 


1.496 x 


10 14 


28 


285.7- 


298.5 


9.076 x 


10 14 


29 


298.5 - 


303.0 


3.217 x 


10 14 


30 


303.0 - 


307.7 


4.484 x 


10 14 


31 


307.7 - 


312.5 


4.857 x 


10 14 
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Table 12— Continued 



Number 


Band 


(nm) 


Flux (photons cm 2 s 1 ) 


32 


312.5 


- 317.5 


5.491 x 10 14 


33 


317.5 


- 322.5 


6.022 x 10 14 


34 


322.5 


- 337.5 


2.286 x 10 15 


35 


337.5 


- 357.5 


3.333 x 10 15 


36 


357.5 


- 377.5 


3.717 x 10 15 


37 


377.5 


- 397.5 


3.947 x 10 15 


38 


397.5 


- 547.5 


6.724 x 10 16 


39 


547.5 


- 735.0 


9.806 x 10 16 



